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ABSTRACT 

In  a  preceding  report  [AFOSE  -  1157,  J'^^   I961]      some 
model  Hamiltonians  were  proposed  for  quantum-mechanical  many-body 
systems  with  pair  forces .   For  infinite  systems  in  thermal  equili- 
brium, they  led  to  temperature-domain  propagator  expansions  which 
were  formally  summable  and  expressible  by  closed  equations.   These 
expansions  were  identical  with  infinite  subclasses  of  terms  from  the 
propagator  expansion  for  the  true  many-body  problem.   The  two 
principal  models  corresponded  to  ring-  and  ladder-diagram  s-ummations 
from  the  true  propagator  expansion,  augmented  by  infinite  classes 
of  self -energy  corrections .   The  model  Hamiltonians  were  called 
stochastic  because  they  contained  parameters  whose  phases  were  fixed 
by  random  choices .   In  the  present  paper,  more  general  models  are 
formulated  which  yield  formally  summable  propagator  expansions  for 
finite  systems.   The  analysis  is  extended  to  correlation  and  Green's 
functions  defined  for  nonequilibrium  ensembles .   The  nonequilibrium 
treatment  is  developed  in  the  Heisenberg  representation  in  such  a 
way  that  unlinked  diagrams  do  not  arise.   A  basic  convergence 
question  associated  with  the  formal  closed  equations  for  the  model 
propagators  and  correlation  functions  is  examined  by  means  of  finite- 
difference  integration  of  the  Heisenberg  equations  of  motion.   This 
procedure  appears  to  converge  independently  of  whether  the  pertur- 
bation expansions  for  the  propagators  and  correlation  functions 
converge.   It  yields  substantial  support  for  the  validity  of  the 
formal  closed  model  equations. 


1 .   IIWRODUCTION 

In  a  preceding  paper   (cited  herein  as  l),  some  so-called 
stochastic-model  Hamiltonians  were  developed  for  fermion  or  boson 
many-body  systems  with  pair  forces .   The  models  yielded  temperature- 
domain  particle-propagators  whose  linked -diagram  expansions  could 
be  formally  summed  and  expressed  by  closed  integral  equations. 
The  model  propagator  expansions  were  identical  with  infinite  classes 
of  terms  from  the  propagator  expansion  for  the  true  many-body 
Hamiltonian.   The  model  Hamiltonians  were  called  stochastic 
because  they  contained  infinite  numbers  of  parameters  whose  phases 
were  fixed  by  random  choices.   Certain  of  the  models,  called  the 
ladder  and  ring  models,  had  important  boundedness  properties  in 
common  with  the  true  many-body  Hamiltonian.   For  this  reason,  it 
was  felt  that  they  might  represent  useful  approximations  to  the 
behavior  of  the  true  system,  particularly  at  low  temperatures. 

The  models  described  in  I  led  to  closed  propagator  equations 
only  in  the  limit  of  an  infinite  system.    In  the  present  paper,  we 
develop  more  general  models  which  yield  formally  summable  propa- 
gator expansions  for  systems  of  any  size.   With  the  new  formulation, 
no  special  considerations  are  required  to  treat  boson  systems  at 
very  low  temperatures,  a  case  which  was  excluded  in  I.  Moreover, 
the  entire  procedure  of  obtaining  the  closed  propagator  equations 
can  be  carried  out  in  a  neater  and  more  satisfactory  fashion.   In 
the  limit  of  an  infinite  system,  the  final  propagator  equations 
obtained  from  the  old  and  new  models  are  identical. 


The  models  employed  in  I  were  constructed  by  altering,  in 
a  stochastic  fashion,  the  true  interactions  among  the  momentum 
modes  of  a  second-quantized  fermion  or  boson  field.   Corresponding 
models  for  distinguishable  particles  were  constructed  by  altering 
the  interaction  among  pairs  of  particles  in  a  way  that  was  different 
for  each  pair.   In  the  present  paper,  we  start  with  a  collection  of 
M  similar  fermion  or  boson  systems,  where  before  we  treated  a  single 
many-body  system.   The  M  systems  are  assumed  to  occupy  the  same  space 
but  not  to  interact.   The  particles  in  each  system  are  distinguishable 
from  those  in  all  the  other  systems.   As  a  preliminary  to  constructing 
stochastic  models,  we  introduce  a  collective  description:  We  define 
M  quantized  fields  which  are  linear  combinations  of  the  quantized 
fields  of  the  M  many-body  systems .   Then  we  alter  in  a  stochastic 
fashion  the  true  dynamical  couplings  among  the  collective  fields. 
As  a  result  of  the  alteration,  the  individual  systems  in  the  collection, 
which  are  independent  in  the  true  problem,  turn  out  to  be  dynamically 
coupled  in  the  models.   The  eventual  closed  equations  for  particle- 
propagators  are  obtained  by  considering  a  grand  canonical  ensemble  of 
collections  and  taking  the  limit  M  -*  «•  . 

We  shall  be  principally  concerned  only  with  fermions  and 
bosons  in  the  present  paper.  For  completeness,  however,  we  shall 
describe  in  Appendix  A  the  corresponding  models  for  distinguishable 
particles. 

The  derivation  of  the  closed  propagator  equations  for  the 
models,  both  in  I  and  in  the  present  paper,  involves  a  fundamental 


convergence  question  concerning  the  contribution  of  classes  of 
infinite-order  diagrams.  No  attempt  to  resolve  this  question  was 
made  in  I .   In  the  second  part  of  the  present  paper,  we  investigate 
the  convergence  question  by  extending  the  analysis  to  systems  not 
in  statistical  equilibrium  and  treating  equilibrium  as  a  limiting 
case.  The  non-equilibrium  formalism  which  we  use  is  a  direct 

adaptation,  to  second-quantized  fields,  of  a  method  previously 

2 
developed  for  the  theory  of  turbulence.    Instead  of  the  temperature- 
domain  particle  propagators,  we  employ  correlation  and  Green's 
functions  which  are  defined  for  time -dependent  statistical 
ensembles  where  no  temperature  exists.   In  the  limit  of  thermal 
equilibrium,  these  quantities  are  related  by  analytic  continuation 
to  the  temperature -domain  propagators . 

As  we  shall  see,  there  exist  both  primitive  and  irreducible 
linked-diagram  expansions  which  give  the  evolution  of  the  correlation 
and  Green's  functions  forward  in  time  from  a  given  initial  instant. 
It  is  convenient  to  develop  the  expansions  in  the  Heisenberg  repre- 
sentation.  Then,  unlinked  diagrams  do  not  arise,  and  therefore 
need  not  be  eliminated.   As  in  the  equilibriim  case,  the  stochastic 
model  Hamiltonians  yield  irreducible  linked-diagram  expansions  which 
are  formally  siramable  and  expressible  by  closed  integral  equations. 

In  the  time -dependent  treatment,  however,  the  validity  of 
the  formally  closed  equations  can  be  examined  by  methods  which  do 
not  involve  perturbation  analysis.   Our  procedure  is  to  replace 
the  exact  Heisenberg  equations  of  motion  by  corresponding  difference 
equations  involving  small  time-increments .   Evidence  is  presented 


that  this  procedure  should  yield  convergent  results  for  the  correlation 
and  Green's  functions,  as  the  increment- size  is  decreased  to  zero, 
independently  of  whether  the  perturbation  expansions  for  these  functions 
converge. 

The  finite-difference  approach  yields  what  we  regard  as 
substantial  support  for  the  validity  of  the  formal  closed  model 
equations.   It  should  be  stressed  at  the  outset,  however,  that  our 
analysis  is  neither  complete  nor  rigorous  and  therefore  is  not  conclusive. 

2.   GENERALIZED  STOCHASTIC  MODELS  FOR  FERMIONS  AND  BOSONS 

2.1.   Collective  Representation  of  a  Collection  of  Systems 

Let  us  consider  a  system  of  fermions  or  bosons  with  Hamiltonian 


^0  "  ^k^k^k  ^k  '  ^^'^^ 

\  =   Kprs\-s\+p,r+s^kW^s  '  ^^'^^ 


for  a  particle  of  momentum  k,  the  e,  are  the  free-particle  energies, 
V  is  a  Fourier  component  of  the  pair  potential,  and  the  summation  is 
over  all  momenta  allowed  by  cyclic  boundary  conditions  on  the  walls  of 
a  box  of  volume  Q  =   L  . ■  We  take  ^  =  1.    This  is  identical  with  the 
Hamiltonian  I:(5.l),  I: (5.2).^  We  shall  continue  to  call  it  the  true 
Hamiltonian.  As  in  I,  the  pair-potential  in  x  space  is  given  by 

V(x)  =  E^Vj^exp(ik-x)  (2.5) 


and  obeys 

V(x)    =  V(-x),        Y^  =  Y_^   ,        V^  =  V^    .  {2.h) 

We  require  fiirther  that  V(x)   be  finite  everywhere  and  that 

V^  ^  0(k"^),  k  -  «    .  (2.5) 

Now,    instead  of  a   single  many-body  system,    let  us   consider 
a  collection  of  M  similar   systems  described  by  the   individual 
second-quantized  fields 

tr^jCx)    =  "'%  q.k[n]    exp(ik-x)  (n  =  1,    2,  ...,M).  (2.6) 

Let  the  true  total  Hamiltonian  for  the  collection  of   systems  be 

H  =  E      Hr  -]    , 
n      [nj    ' 

^[n]    =  \  \^k[n]\[n]  (2.?) 

+  22:kprs\-s\+p,  T+s\  [n]  \  [n]  \  [n]  ^^s  [n]  ; 

We   shall  take  the   commutation  rules  as 

W[n\>    S[m]^    +  =   °'  ^^k[nj    >    S  M '^    +  ^   ''"^^  *    ^^'^^ 

Thus  the  individual  systems  exist  in  the  same  box  but  othervfise  are 
entirely  independent.    In  (2.8),  and  throughout  the  paper,  the 
upper  sign  of  a  double  sign  refers  to  fermions,  and  the  lower  to 
bosons . 

As  a  preliminary  to  constructing  stochastic  models,  we 
shall  express  H  in  terms  of  collective  fields  which  are  linear 


combinations  of  the  \(/r  •](x).   Let  us  restrict  M  to  the  form 
M  =  2S  +  1,  where  S  is  a  positive  integer.   We  define  the 
collective  fields  ij/  (x)  by 


tjx)  =  m'^Z^  exp(i2jran/M)i^r  i(x)  , 
tj(x)  =  m"2z^  exp(-i2rtan/M)^r  i+(x) 


(2.9) 


where 

a  =   -S,  .  .  .,-1,0,1. . . .,S    . 
¥e   shall  call  a  a   collective   index.        By  using  the   identities 

m'^Z^  exp[i2rt(a  -   p)n/M]    =   5^p    ,  (2.10) 

m""^Z     exp[i2rta(n  -  m)/M]    =  6        , 
we  may  invert    (2,9)   and  obtain 


tr  -|(x)   =  m"2z^  exp(-i2rtan/M)\l/^(x)    .  (2.11) 


q,    "*"  may  be  defined  by 
ka 


tjx)    =  n"2Zj^   q^^  exp(ik-x)    ,  (2.12) 


\a=  ^'"^n  exp(i2mn/M)qj^^^^ 


'lk[n]    =  ^'"^a  exp(-i2rtan/M)qj^^ 


(2.15) 


It  is  clear  from  this  that  the  transformation  to  the  collective 
fields  represents  a  complex  rotation  in  the  space  of  the  q  r  -i  . 
By  (2.8),  (2.10),  and  (2.13),  the  commutation  rules  in  the 
collective  representation  are 

Let  us  adopt  the  cyclic  convention 

a=a+M,  n=n+M  (2.15) 

for  collective  indices  and  individual-system  labels.  This  clearly 
is  consistent  with  ( 2.9)-(2.1l) .   Then,  by  using  (2.13)  and  (2.10). 
we  may  rewrite  H  in  the  form 


So  =  "a'^k^k'ika^ka  ' 


(2.16) 


h  =  *^'\p^iX^kprs\-s\+p,r+sVp,^+xV\3V^^s;.  ,  (2-1^) 


where  5        is  to  be  interpreted  by  the  cyclic  convention.   We 
see  from  (2.17)  that  the  collective  fields,  in  contrast  to  the 
individual- system  fields,  are  dynamically  coupled  by  the  interaction 
V(x). 

It  is  clear  from  the  definitions,  and  also  from  an  exami- 
nation of  (2.1^+)  and  (2.17),  that  the  momenta  k  and  the  collective 
indices  a   have  some  formal  properties  in  common.   In  particular. 
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the  factor  5  „   ,  in  (2.17)  gives  a  conservation  property 

CL+p  ,  l^+A 

analogous  to  momentum  conservation.   The  imderlying  similarity  may 
be  expressed  as  follows.    Consider  the  case  where  Q.   is  very  large, 
so  that  each  many-body  system,  described  by  a  field  ^r„"](x)'  i^ 
effectively  composed  of  many  spatially  localized  sub-systems.   The 
the  Fourier  components  q   ,  considered  as  functions  of  k,  provide 
a  collective  description  of  the  sub-systems  in  x  space.   On  the 
other  hand,  suppose  we  take  M  large  instead  of  Q.   large.  Then 
the  q.,  ,  considered  as  functions  of  a  ,  provide  a  collective 
description  in  "n  space." 
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2.2.  Formulation  of  Models 

In  order  to  construct  stochastic  models,  we  replace  (2.17) 
by  the  more  general  form 

2i  =  *^"\pnX^prs\-s^aPn). 

^  \+p,r+s^a+3,n+X   '^ka'^Sp^V'^sX-  (2.l8) 


The  c  number  coefficients  $  „  ^^  play  a  role  analogous  to  that  of 
the  $.     in  I:{j>.h).      We  impose  upon  these  coefficients  the 
Hermiticity  and  symmetry  constraints 

apuX  =  -A|jpa*  '  ^apiaX  =  ^paA|i  '  (2-19) 


Then  we  make  stochastic  assignments  of  values  to  the  $  o  >  so  as 
to  obtain  the  several  stochastic  models.   We  give  below  the  pre- 
scriptions which  yield  the  present  analogs  of  the  ladder,  ring, 
random-coupling,  and  Hartree-Fock  models  defined  and  discussed  in  I. 

In  each  model,  the  form  of  $  „  ,  is  identical  with  that  of  $, 

aPi-iA  kprs 

in  I: Sec.    J. 

Ladder  Model 


'ap^X  "   — ±^L-v   "p^    ■    -^^,j    ,  (2.20) 


'     a         Qq  (2.21) 

ap  =     Pa   .  ^  ■' 


Determine  the  real  phases  6   by  independent  random  choices  in  the 
interval  (O,  2n) ,    subject  only  to  (2.21). 
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Ring  Model 
Take 

%^X  =   ^^P[i(V^  ^(B^)J'  (2.22) 

^oA  =  -  V  •  (2.23) 

Determine  the  phases  6  i^by  independent  random  choices  subject  only 
to  (2.25). 

Random- Coupling  Model 
Take 

®ap|aX  "  '^Xupa  '   ®aPiaX  "  ®paA|a  '   ®aP|aX  "  ®pa^X  "      (2-25) 
Determine  the  phases  6  „  ,  by  independent  random  choices  subject  only 


to  (2.25). 


Ifertree-Fock  Model 
Take 


^  CO      =  <^  o  o  =  1  '     "J"  a  ^  =  0    (a  +  u  or  X)  .      (2.26) 
appa    apap  aPiaX        v    h-     /       \  / 


There  are  no  random  parameters . 

The  models  obtained  by  the  prescriptions  just  given  are 
stochastic  in  the  same  sense  as  discussed  in  I: Sec. 2. 2.   The  values 
of  the  f  „  ■,  are  determined  by  random  choices.   Once  obtained,  however, 
the  values  are  fixed,  and  we  work  thereafter  with  the  definite  Hamil- 
tonian  embodying  these  values.   To  obtain  an  insight  as  to  the  physi- 
cal interpretation  of  the  present  models,  it  is  of  value  to  transform 
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(2.18)  back  to  the  individual-system  representation.   We  find 

H.  =  i-Z    ,  ,E,    V,   Ar   ,  ,-| 
—     i„,T  ]jpj,3  j^_g  ^nmm'n'J 

^  \+p,r+s\[n]\[m]\[m-]S[n']    '      ^^-S?) 


nmn ' 


where 


A  r  -1    =  m"  Z  8  $ 

[nmm'n'J  aPiaX     a+p,n+A.  aPiaX 


X  exp[i2rt(-na-mp-m'p.+n'>L)/M]    .  (2.28) 


When  all  the   ^^   -^  =  ^^    i't  follows  from   (2.10)   that 


l_nmm  n  'J 


Then  H  reduces  to  (2.7),  giving  us  back  the  true  problem  in  which 
the  individual  systems  are  dynamically  independent.   In  the  models, 
however,  Ar   ,  ,-|  is  not  diagonal,  and  the  systems  in  the  collec- 
tion are  coupled.   They  exchange  particles  as  well  as  energy  and 
momentum.   We  see  from  (2.27)  that  the  elementary  interaction  is  a 
collision  in  which  a  pair  of  particles  from  systems  n'  and  m' 
are  destroyed  and  a  pair  in  systems  n  and  m  created. 

2.3  Bounds  on  the  Ladder  and  Ring  Model  Hamiltonians 

Bounds  on  the  eigenvalues  of  H  for  the  ladder  and  ring  models 
may  be  derived  in  close  analogy  to  I: Sec. 3.   In  the  case  of  the  ladder 
model,  let  us  define  the  effective  two-body  amplitudes 

Xjx',x)  =M-*Z^5^^^^^yx')t^(x)exp(ie^^)  ,  (2.29) 


Ik 


where  (2.I5)  is  to  be  used  in  interpreting  5  ,   .  Then,  noting 
(2.12),  we  may  write  (2.l8)  for  the  ladder  model  in  the  form 


-i  =i^e  /7v(x-x')X^+(x',x)Xjx',x)d^xd\'  .     (2.50) 

It  is  clear  from  (2.30)  that  H.  is  a  non-negative  operator  if 
V(x)  =  |v(x) I  everywhere.   It  then  follows  that  the  eigenvalues 
of  H  are  all  non-negative  in  this  case. 

In  the  case  of  the  ring  model,  we  define  the  effective  density 
components  p   hy 

P^(x)=Zj^Pj^^exp(ik-x),  (2.31) 

p^(x)  =  M-\V(,)^^^^(,)exp(i^^^^^P  .  (2.32) 


Then,  using  (2.1*+)  and  (2.23),  we  find 

where  N,  the  operator  for  the  total  number  of  particles  in  the 
collection,  is  defined  by 

—    a  a    1^   L^J 
\=^%a\a^    ^[n]  =  ^k  ^^k  [nj   ^k  [n]  ' 

If  V,  =  |V  I  for  all  k,  then  the  first  term  on  the  right  side 

of  (2.33)  is  a  non-negative  operator,  and  the  energy  per  particle 

in  any  eigenstate  of  W  is  bounded  from  below  by  -2^(0). 
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5.   PROPAGATOR  EXPANSIONS  FOR  THE  MODELS 

In  analogy  to  I:Sec.i|,  let  us  define  the  temperature-domain 
propagator  \rnm-|(u,u')  by 

^kH^-'-')  =  -  <'\[n]^^K[raf(-">y       '  (5.1) 

where,  for  any  operator  B, 

B(u)  =  exp(uH)Bexp(-uH),    B  (u)  =  exp(uH)B  exp(-uH)  .   (3.2) 

The  ordering  operator  T  is  the  same  as  in  I,  and  the  brackets 
denote  an  average  over  a  grand  canonical  ensemble  of  collections 
of  systems:   For  any  B, 

<B>=  Tr  {exp[-p(H-nN)]B}  /Tr^xp[-p(H  -  uN)]}  .      (3.5) 

Let  us  also  define  the  propagator  in  the  collective  representation 
by 

\ap^^'-'^  =  -KJ-Kp'(-')>    •  (5.M 

'Energy '-domain  propagators  S^^  r^J  (  ^g^ )  and  S   (^^)  may  be  defined 
in  terms  of  Sj^  rj^l  (u,u  ' )  and  Sj^(^n(u,u')  by  means  of  relations  of 
the  form  I: (4.6) . 

For  free  particles  (H^  =  O),  it  follows  from  (2.8),  (2.1^)-), 
and  (2.16)  that 

S,  r  -i(M  =  5  S,  (°Vf  )  , 
kLnmJ  ^^a'    nm  k   ^^a'^  ' 

,  ,  (3.5) 
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where  S^    '(^    )    is  given  by  I: (4.11).        For  the  true  Hamiltonian 
(^^^^ap^X  =  D'   ^e  have 


R  r    n(^  )  =  5    a  (^  )  ,  (3.6) 

k  LnmJ  ^  ^a  -^  ran  k  ^  ^a ' 

where  S  (^  )  is  the  true  propagator,  as  defined  in  I,  for  any  one 
of  the  systems  in  the  collection.  This  follows  immediately  from 
the  independence  of  the  individual  systems.   In  this  case,  we  find 


where  we  have  used  (2.15)  and  (2.10). 

In  correspondence  to  I: (4. 8)  and  I: (4.10),  the  mean  number 
of  particles  and  the  mean  energy  in  the  collection  may  be  expressed 
for  the  models  in  the  forms 


N(P,^)  =  ±^~\^^\^iQ--v{^,^) 


(5.8) 


(3.9) 


=  tP''^nkaS^[nn](^a)-P(^a^) 


E(P,^)  =  l*P''z^JV^a)SkcKx(^a)-P(^a^) 

=  t^'\ka(V^a)\[nn](^a^-P(^a^) 


Equation  (3-9)  may  be  derived  in  the  same  way  as  I: (4.10).  The 

mean  number  of  particles  per  system  and  the  mean  energy  per  system 
are  then  given  by 

N(p,n)  =  M"^N(P,h),     E(P,h)  =  M"^E(p,n)  .  (3.10) 
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For  the  true  problem  (all  $,     =  1  in  I  and  all  $  „  ,  =  1  in  the 
Kprs  apioA 

present  case),  N(3,|a)  and  E(f3,u)  are  exactly  the  same  functions  of 
3  and  (j  here  as  in  I. 

Let  us  now  consider  the  primitive  linked-diagram  expansion 
for  the  model  S^  w^^a^'   ^i^^t  of  all,  it  is  easy  to  verify  that 

Skap(^a)  =  Vkoa^^a)  '  (5.11) 

This  follows  immediately  from  (3-5)  and  the  presence  of  the  conser- 
vation factor  5      ,  in  (2.l8).   Because  of  the  similarity  of  the 
ways  the  momenta  k  and  the  collective  indices  a  enter  H,  the  primitive 
linked-diagram  expansion  for  S   (^  )  is  given  by  rules  which  are 
obvious  generalizations  of  rules  1-8  of  I: Sec.  k.     We  shall  state 
here  only  the  changes  which  are  required  in  the  rules  of  I: 


Rule  5.   Label  the  external  lines  with  collective  index  a  as  well  as 
momentum  k.   Label  the  internal  lines  with  collective  indices  a',  a" 
as  well  as  with  momenta  k',  k  ",  .  .  .  . 

Rule  5.  With  each  vertex  labeled  as  in  Fig.  1  of  the  present  paper, 
associate  a  factor 

-m"   6~  $  5  V        5  6 

P7IJ.X   p+7,M.+X  p-s  p+q, r+s  b+c,d+e 

Riile  7-      Sum  over  all  the   intermediate  collective  indices  a',   a"  ,... 
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as  well  as  over  all  the  intermediate  momenta  and  'energies, 


According  to  the  present  rules,  the  contribution  of  any  diagram 

to  S,   (?)  has  the  following  form.   It  consists  of  a  product  of  V 
koa  a ' 

and  S^   factors,  summed  over  momenta  and  energies,  which  is  multi- 
plied by  a  product  of  $ 's,  summed  over  collective  indices.   The 
latter  summed  product  contains  the  entire  dependence  of  the  contri- 
bution upon  the  ^'s   and  upon  a.   The  summed  product  of  V  and  S 
factors  is  identical  with  the  contribution  which  the  diagram  makes 
to  S  (?  )  in  the  true  problem,  where  S  (?  )  is  defined  as  in  I. 
In  order  to  express  this  result  compactly,  let  us  write  S,  (C  )  for 
the  true  problem  in  the  form 

S,  (?  )  =  a  ^°^?  )   +  ^"  ,^     r   (k,a)  .         (3.12) 
k^^a'^    k   ^^a-^    n=l  p  n;p^  '  ' 

Here  T   (k.a)  is  the  contribution  from  a  particular  distinct 
n;p 

linked  diagram  with  n  vertices.   The  index  p  takes  the  Aralues 
1,  2, ...,R(n),  where  R(n)  is  the  number  of  distinct  linked  diagrams 
with  n  vertices;  it  labels  the  R(n)  diagrams  according  to  any  con- 
venient scheme.   We  may  now  write  the  primitive  linked-diagram 
expansion  for  the  model  propagator  S  ^^(^g^)  ^.s 

(0), 


s, 


(?  )  =  s/^^(0  +  S  ,E   C    (a)r   (k,a)  ,    (3.13) 
L^^a'^    k   ^^a'    n=l  p  n;p^   n;p^ 


where  C   fa)  denotes  the  summed  $  product  associated  with  the  pth 
distinct  diagram  with  n  vertices.   The  T^  (k,a)  in  (3-13)  are 
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precisely  the  same  quantities  as  in  (5.12). 

We  may  illustrate  the  results  stated  in  the  last  paragraph 
by  considering  the  contributions  to  S,   (^  )  made  by  two  simple 
primitive  diagrams .   Let  us  assign  the  value  p  =  1  to  the  diagram 
of  Fig.  2a.   Then  the  rules  for  constructing  diagram  contributions 
give 

r^.^(l.,a)  -   +3%  V„S^(°'(C„)s/°'(C,)exp(t,6)sJ°>(C^)  ,    (^.A) 

C^.^(a)  =  M-^Ep*_^^^  ^  (3.15) 

Let  us  assign  the  value  p  =  2  to  the  diagram  of  Fig.  2b.   Then  the 
rules  give 


r^^^(.,a)  =  f\^^^   \-s\-p\'°*«a)S'°'«b> 


X  ^'°'(^.)^..p-.""(^a*-c)\'"'(^a>  '      C-*^) 


_2 


Equations  (3.15)  and  (3.1?)  illustrate  a  general  property  of  the 
C  ,  (a).   The  niunber  of  M   factors  is  always  equal  to  the  number 
of  free  indices  in  the  summation.   Hence  in  the  true  problem  (all 
^       ,  =  ]),  we  have  C   (a)  =  1  for  all  n  and  p,  as  we  must  for 
consistency. 

Let  us  now  consider  the  limit  M -*  °°,  which  represents  an  infinite 
collection  of  many-body  systems.   We  keep  (i  fixed  as  we  take  the  limit. 
For  the  true  problem,  N(P,h)  and  Y,{^,\i-)   are  then  independent  of  M, 
as  we  have  noted  previously.   It  will  appear  from  what  follows  that 
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W(P,|a)  and  E(|3,h)  for  the  models  are  independent  of  M  in  the  limit. 
The  present  limit  M  ^  oo  plays  the  same  role  as  the  limit  fi  -*  <» 
(infinite  volume)  did  in  I .   We  foimd  the  following  general  result 
in  I:   For  any  given  model,  a  given  diagram  survived  in  the  limit 
Q  ->■  00  if  and  only  if  the  associated  $  product  had  the  value  one  for 
all  values  of  the  intermediate  momenta  permitted  by  momentum  con- 
servation.  The  corresponding  result  in  the  present  case  is  that  a 
given  diagram  survives  in  the  limit  M  ->■  «>  if  and  only  if  the  associated 
$  product  [i.e.,  the  summand  of  C  ,  (a)]  has  the  value  one  for  all 
values  of  the  intermediate  collective  indices  permitted  by  collective- 
index  conservation.   For  each  surviving  diagram  we  therefore  have 
C  _  (a)  =  1,  as  in  the  true  problem.   The  demonstration  of  the 
present  result  is  neater  than  that  in  I  because  the  entire  collective- 
index  dependence  of  each  diagram  contribution  is  contained  in  the 
$  product,  whereas  in  the  former  treatment  the  <J>  factors^  S    factors, 
and  V  factors  all  depended  on  the  intermediate  momenta. 

The  correspondence  between  the  survival  of  diagrams  in  the 
old  treatment  and  the  present  should  become  clear  upon  comparing 
the  following  two  examples  with  the  corresponding  examples  in 
I: Sec. 5. 3.   For  the  ladder  model,  we  have 

=  1,  (5.18) 

where  we  have  used  (2.20)  and  (2.21).   Thus  this  diagram  survives 
and  gives  the  same  contribution  as  it  does  in  the  true  problem.  For 
the  ring  model,  however,  we  have 


c^J^)  =  M"'Sx^^Pti(e^+%,^^p_^)  +  i(e,3+e^^p_,^J]  ,    (5.19) 


by  (2.22).   Except  for  special  values  of  the  indices  3  and  X, 
the  phase  of  the  summand  does  not  vanish  as  a  result  of  (2.25),   Con- 
sequently, for  a  typical  assignment  of  the  phases  6   by  random 
choices,  we  have  Cp  „(a)  -»  0  as  M  -*  <»,  and  the  diagram  does  not  sur- 
vive.  The  detailed  reasoning  follows  that  in  I. 

Similar  analysis  may  be  applied  to  any  diagram  of  any  finite  order. 
It  may  thereby  be  verified  that  the  surviving  diagrams  for  each  of  our 
present  models,  in  the  limit  M  -»  «,  are  precisely  the  same  as  those 
for  the  model  of  corresponding  name  in  I: Sec. 5,  in  the  limit  fl  -♦  oo. 
This  is  a  consequence  of  the  fact  that  the  rules  for  diagram- contributions 
in  the  two  cases  yield  $  products  with  the  same  structure. 

We  shall  now  make  a  fundamental  ass-umption  which  corresponds 
to  assumption  5  of  I:Sec.  5'1»  We  assume  that  diagram- classes 
which  do  not  sui^ive  in  any  finite  order  make  no  contribution  to 
S,   (t  ),  in  the  limit  M  ^  oo,  when  they  are  summed  to  all  orders. 
The  implications  of  this  assumption  will  be  discussed  in  Sec.  h, 
and  in  Sec.  7  we  shall  outline  a  partial  justification. 

Since  the  diagrams  we  have  identified  as  surviving  all  have 
C  _  (a)  =  1,  an  immediate  consequence  of  our  assumption  is  that  S,   (^  ) 
is  independent  of  a  in  the  limit.   Thus  we  may  define  S   (^  )  for 
the  models  by 

\^U=\aa(^a)        ^^  ^  ^^  '  ^5.20) 

Then,  by  (2.13),  (5.11),  and  (2.10),  we  have 
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S,  r  t(C  )  =  5  a  (^  )        (M  -  ")  .  (3.21) 

In  the  true  problem,  these  relations  hold  for  any  M. 

Let  us  now  consider  the  irreducible  linked- diagram  expansion 
for  S  (^  ),    as  defined  by  (3.20),  and  for  R  (^  ),  as  defined  in  terms 
of  this  S,  (^  )  by  I: (4.12).   These  expansions  are  obtained  by  amend- 
ing the  rules  for  the  primitive  diagram  expansion  vith  rules  2',  h' , 

4"  ,  and  8'  of  I: Sec. 4.   Since  each  S   (^  )  is  independent  of  a  in 

koxx  a 

the  limit  M  -*  oo,  it  follows  that  the  entire  collective-index  dependence 
of  any  diagram- contribution  in  the  irreducible  expansions  is  contained 
in  the  associated  0  product  and  therefore  is  given  by  the  same  quan- 
tity C  .  (a)  associated  with  that  diagram  in  the  primitive  expansion. 
From  this  it  follows  that  the  surArLval  rules  for  diagrams  in  the 
irreducible  expansions  are  precisely  those  for  the  primitive  expansion: 

A  diagram  survives  if  and  only  if  the  summand  of  C   (a)  has  the  value 

n;p 

one  for  all  values  of  the  intermediate  collective  indices  permitted 
by  collective-index  conservation.   This  same  conclusion  may  also  be 
reached  by  reasoning  of  the  kind  employed  in  I: Sees.  5.2  and  5«5 
concerning  recovery  of  the  primitive  expansion  from  the  irreducible 
expansion. 

The  results  stated  in  the  preceding  paragraphs  may  be  combined 
to  give  the  following  conclusion:   In  the  limit  M  -*  «>,  only  those 
classes  of  terms  surArLve,  in  the  irreducible  expansions  for  S,  (^^) 
or  H  (^  ),  which  were  identified  as  surviving  for  the  corresponding 
models  in  I: Sec.  5.   In  each  model,  the  surviving  terms  are  identical 
in  form  with  the  corresponding  tenas  in  the  irreducible  expansion  for 
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the  true  problem.   Thus,  for  each  of  our  present  models,  we  are 
led  to  precisely  the  same  final  closed  equations  for  S  (^  )  as  we 
obtained  in  I: Sec.  5.   The  essential  difference  is  that  now  these 
equations  are  obtained  for  finite  Q.      Asstmiptions  5  and  k   of  I:Sec.  5.I 
have  not  been  made  here.  An  important  consequence  is  that  our  present 
resiolts  apply  to  boson  systems  at  low  temperatures.   A  further  advan- 
tage of  the  present  procedure  is  that  the  double  limiting  process 
employed  in  I  for  long-range  potentials  is  now  unnecessary. 

h.      CONVERGENCE  PROPERTIES  OF  THE  MODEL  PROPAGATOR  EXPANSIONS 

For  each  stochastic  model,  let  us  diAd.de  all  the  primitive 

linked  diagrams  into  two  classes,  wanted  and  -unwanted.  Wanted 

diagrams  are  all  those  which  we  identified  in  Sec.  5  as  surviving 

in  the  limit  M  -*  «.   The  $  products  associated  with  them  have  the 

value  one  for  all  values  of  the  intermediate  collective  indices 

permitted  by  collective -index  conservation.   Consequently,  we  have 

C  _  (a)  =  1  for  these  diagrams.   The  unwanted  diagrams  are  all  the 

others.    The  $  products  associated  with  them  have  stochastic  phases 

for  all  but  special  values  of  the  Intermediate  collective  indices. 

Consequently,  we  have  C   (a)  -*  0  (M  ^  ~)  for  any  given  unwanted 
n;p 

diagram.   We  shall  examine  very  soon  the  strength  with  which  the 
unwanted  C   (a)  vanish  in  the  limit. 

The  fundamental  ass-umption  made  in  Sec.  5  concerning  the  con- 
tribution from  all  unwanted  diagrams  may  be  expressed  in  the  form 

l"   ^LC        (a)  r   (k,a)  ->  0     (M  ^  «)  ,  (^.l) 

n=l  p  n;p^  '   n;p^ 
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where  E  denotes  a  summation  over  imwaxited.  diagrams  only.   The 
essential  point  in  (^.l)  is  that  we  must  take  the  limit  n  -»  «  before 
the  limit  M  ^  oo.   If  we  did  not,  we  woijld  not  be  employing  the 
full,  formally  exact  perturbation  series  for  the  model,  and,  con- 
sequently, we  would  be  imable  to  assert  that  the  results  embodied 
the  boundedness  properties  of  the  model  Harailtonian.  We  wish  to 
make  it  clear  in  the  present  Section  that  the  validity  of  (^.1)  is 
far  from  obvious.   We  shall,  in  fact,  obtain  the  negative  result 

Cl^p  l%p(-)^n:p(^'^)l  =  ~       (^-^  ^)  •         ^^-2) 

This  result  does  not  mean  that  (^.1)  cannot  be  true.   However,  it 

shows  that  (^.l)  can  hold  only  if  there  are  cancellations  among 

the  contributions  from  different  diagrams.   We  shall  not  attempt 

to  investigate  the  latter  question  directly.   Instead,  we  shall 

present  in  Sec.  7  some  evidence  for  the  validity  of  (^.1)  which  is 

not  based  on  perturbation  theory. 

In  order  to  establish  (^.2),  we  must  estimate  three  things: 

the  strength  with  which  the  unwanted  C   (a)  vanish  for  large  M, 

n;p 

the  number  R'(n)    of  distinct  unwanted  diagrams   for  large  n,    and  the 

magnitude  of  the  true-problem  diagram  contributions  r   .    (k,a)    for 

n;p 

large  n.   Let  us  consider  the  C  ,  (a)  first.   In  every  C  ,  (a), 

the  number  of  independent,  summed  indices  given  by  index-conservation 

and  the  number  of  M~   factors  are  both  equal  to  n.    The  $  product 

which  comprises  the  summand  of  any  C  ,  (a)  has  unit  modulus  for 

all  of  our  models  but  the  Hartree-Fock.    For  that  model  the  0  product 

is  either  zero  or  one.   It  follows  that,  for  all  the  models,  the 
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the  C  .  (a)  all  satisfy 

There  are  two  types  of  contributions  to  the  unwanted  C  .  (a). 
First,  in  each  of  the  models,  there  are  certain  restricted  values 
of  the  intermediate  collective  indices  for  which  the  0  product  is 
identically  one.   (Cf.  the  discussion  in  I:Sec.  5.2.)  We  may  estimate 


Hence  the  term  in  the  sum  with  all  collective  indices  =  a  has  the  value 
M~  ,  which  then  constitutes  a  lower  bound  to  the  magnitude  of  the 
contribution.  An  upper  bound  may  be  given  in  the  form  KM  ,  corres- 
ponding to  the  restriction  of  just  one  of  the  intermediate  indices 
to  K  special  values. 

The  second  type  of  contribution  comes  from  the  terms  in  the 
simi  which  have  stochastic  phases.   It  arises  in  all  the  models  but 
the  Hartree-Fock.    There  are  M^  terms,  each  of  unit  modulus,  in 
the  sum  which  (when  multiplied  by  M  )  comprises  C   (a).   If  the 
phase  of  each  tenn  were  fixed  by  an  independent  random  choice,  then, 
for  a  typical  set  of  choices,  the  result  would  be 

C^.  (a)  =  M""o[V^]  =  oCm"""/^)  . 

In  general,  however,  the  JVi  phases  are  not  completely  independent. 

2 
The  independent  choices  are  of  the  0(M  )  linearly  independent  6  „ 

or  e   (ladder  or  ring  model)  or  the  0(M  )  linearly  independent 

e     (random- coupling  model).   In  consequence,  the  contribution 

to  C   (a)  from  the  stochastic-phase  terms  may  vanish  either  more 
n;p 
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or  less  strongly  than  M    .   Furthermore,  if  the  unwanted  primitive 

diagram  contains  a  wanted  diagram  as  a  self -energy  part,  variation 

of  the  intermediate  indices  corresponding  to  this  self -energy  part 

will  not  give  any  change  in  the  phase  of  the  $  product.   This  will 

result  in  C   (a)  vanishing  less  strongly. 
n;p 

On  the  basis  of  the  preceding  paragraphs,  we  shall  assume 
the  following  asymptotic  behavior  for  the  unwanted  C  _  (a): 

C   (a)  -  0,        M-  oo  ,  (i|.4) 


|C^.p(a)l  >M-^      M--  .  (i|.5) 

Two  restrictions  on  the  validity  of  these  relations  should  be  noted. 
We  assume  a  typical  assignment  of  the  independent  random  phase- 
parameters  in  asserting  (^.^).   (Cf.  the  discussion  in  I:Sec.  2.2.) 
In  obtaining  the  lower  bound  (^.5),  we  have  ignored  the  contribution 

to  C   (a)  from  the  terms  with  stochastic  phases.   However,  there 
n;p^ 

may  be  some  special  combinations  of  values  for  n,  p,  and  cc  such 

that  the  non-stochastic  contribution  to  C  .  (a)  is  very  nearly 

n^p 

cancelled  by  the  stochastic-phase  contribution,  with  the  resiolt  that 
(4.5)  is  violated.   For  a  typical  assignment  of  phase-parameters, 
these  cases  constitute  a  set  of  relative  measure  zero  in  the  limit 
M  -»  00.   We  shall  ignore  them  in  discussing  (4.2). 

In  order  to  boimd  R'(n),  the  number  of  unwanted,  distinct 
linked  diagrams  with  n  vertices,  we  consider  a  special  class  of 
linked  diagrams,  which  we  shall  call  linear  diagrams.   Let  us  form 
an  Incomplete  diagram  with  n  vertices  as  illustrated  in  Fig.  3  for 
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the  case  n  =  4.   To  form  the  linear  diagrams,  ve  now  connect 
the  free  outgoing  lines  with  the  free  incoming  lines  in  all  possible 
ways.   There  are  n!  ways  of  doing  this,  and  each  gives  a  distinct 
linear  diagram.   Now  we  note  that  the  linear  diagrams  are  a  very 
restricted  class  of  all  the  distinct  linked  diagrams.   Hence,  for 
large  n,  we  have  R(n)»  n!.  Moreover,  we  note  that  the  unwanted 
linked  diagrams  constitute  most  of  the  linked  diagrams  for  large  n, 
no  matter  which  of  our  models  we  take.   The  wanted  diagrams  in 
each  case  are  very  special  classes.    Consequently,  we  may  safely 
assume  the  bound 

R'(n)   >  n!       (n  ^  «>  )  .  (k.G) 

We  have  finally  to  estimate  the  r  _  (k,a)  for  large  n. 

On  the  basis  of  (2.5),  I: (4.11),  and  I: (4.7),  we  assume  that  for 

any  given  n  and  p  the  summations  over  intermediate  momenta  and 

intermediate  'energies'  converge  at  infinity  (cf.  assumption  2  of 

I: Sec.  5.1).   It  follows  that  r  .  (k,a)  can  be  bounded  in  the  fashion 

n^p 

|b(k,a)|"<  lr^.p(k,a)|  <  |B(k,a)|"  (n  -  «)  ,   (4.7) 

where  b(k,a)  and  B(k,a)  are  parameters  independent  of  n  and  p.   In 
general,  we  will  have  |b(k,a)|  >  0.    We  shall  not  attempt  to  prove 
(4.7)  here. 

Equation  (4.2)  follows  immediately  from  (4.5),  (4.6),  and 
(4.7).   We  have,  furthermore, 

S"  ,E  I  r   (k,a)|  =  00  .  (4.8) 

n=l  p'   n:p^  ■'  '  ' 


28 


That  is  to  say^,  the  primitive  linked-diagram  expansion  for  the  true 
problem  is  not  absolutely  convergent  when  taken  diagram-by -diagram. 
This   result  does  not  preclude  the  weaker  property  that 

E°°   Jz  r        (k,a)| 
n=l'    p  n;p^    '    '  ^ 

converge.   The  latter  would  correspond  to  absolute  convergence  of 
S  (^  )  as  a  power-series  in  X  if  the  interaction  potential  V(x) 
were  replaced  by  XV(x). 

We  have  already  noted  that  (^.2)  does  not  mean  (^.l)  cannot 
be  true.   However,  it  certainly  indicates  that  the  assumption  expressed 
by  (^.l)  must  be  regarded  with  suspicion  in  the  absence  of  supporting 
evidence.  We  shall  attempt  to  develop  such  evidence  in  the  remainder 
of  this  paper  by  turning  to  the  more  general  problem  of  linked-diagram 
expansions  for  non-equilibriiim  statistical  ensembles  which  evolve 
in  time. 
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5.   NON-EQUILIBRIUM  CORRELATION  AND  GREEN'S  FUNCTIONS 
5.1  Summary  of  the  Method 

The   formalism  we   shall  outline  in  Sees.    5  and  6  is   an  ex- 
tension,   to  second- quantized  fields,    of  a  non-equilibrium  method 

previously  applied  to  a  classical  field  problem,  the  evolution  of  the 

2 
correlation  tensor  in  turbiilence  dynamics.   The  changes  required 

in  the  present  case  are  due  to  the  different  degree  of  nonlinearity 
of  the  equations  of  motion  and  to  the  q-number  nature  of  the  fields. 
Neither  characteristic  necessitates  a  drastic  modification.  We  shall 
develop  the  formalism  here  ab  initio,  but  not  with  full  proofs  or  in 
complete  detail.   Oxir  interest  here  is  not  in  non- equilibrium  as  such 
but  in  elucidation  of  the  convergence  question  which  we  stated, 
for  equilibrium,  in  Sec.  k. 

The  non-equilibrium  treatment  will  be  carried  out  in  the 
Heisenberg  representation.   Our  procedure  may  be  summarized  as 
follows.   We  define  Heisenberg  creation  and  destruction  operators 
which  coincide  with  the  Schrodinger  operators  at  an  initial  time 
t  ,  and  we  specify  an  initial  statistical  ensemble  by  a  weighting 
operator  which  is  Gaussian  in  the  Schrodinger  operators.   Then  we 
define  correlation  functions  of  the  Heisenberg  operators,  and  also 
Green's  fimctions  which  give  the  average  response  of  these  operators 
to  infinitesimal  external  perturbations.   The  correlation  and 
Green's  functions  are  expressed  in  terms  of  the  Schro'dinger  opera- 
tors by  iterative  solution  of  the  Heisenberg  equations  of  motion. 
The  results  are  evaluated  by  a  statistical  form  of  Wick's  theorem 
which  is  valid  for  any  Gaussian  weighting  operator.   This  yields 
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linked -diagram  expansions  for  the  correlation  and  Green's  functions. 
An  advantage  of  working  in  the  Heisenberg  representation  is  that 
Tjnlinked  diagrams  do  not  arise  and  therefore  do  not  have  to  be 
eliminated.   Both  primitive  and  irreducible  linked-diagram  expansions 
are  obtained.   The  latter  yield  formally  closed  equations  which  give 
the  evolution  in  time  of  the  correlation  and  Green's  functions  for 
each  of  our  stochastic  models. 

5.2  Quantized  Nonlinear  Oscillator 

In  order  to  minimize  notational  complications^  we  shall 
introduce  the  non-equilibrium  treatment  in  tenns  of  a  simple 
example  which  yields  diagram  expansions  identical  in  structure 
with  those  for  the  many-body  problem.   Consider  the  true  Hamiltonian 

^0=  ^^n  ^[n]"^^[n]  ' 

(5.1) 
'n^[n]'^[nj'^[nj^[nj  ' 

where  V  and  e  are  real  parameters  and  the  commutation  relations  are 

[^[n]^  ^[m]]  +  =  0^     ^\r.Y    ^[m]'l  +  =  ^nm  '     ^5-2) 

The  index  n  =  1,  2,...,   W.   specifies  individual  systems  in  a 
collection,  as  before.   This  Hamiltonian  represents  a  collection 
of  trivially  soluble  quantized  nonlinear  oscillators.   In  the 
fermion  case,  H.  actually  vanishes  identically,  but  this  will  not 
affect  the  usefulness  of  (5.I)  for  our  purposes.    In  terms  of  the 
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diagram  expansions,  the  vanishing  of  H.  is  expressed  by  the  exact 
cancellation  of  the  contribution  of  any  non-exchange  diagram  by 
that  of  an  exchange  diagram;  the  formal  structure  of  the  expansior 
is  unaffected. 

The  general  model  Hamlltoniaxi  corresponding  to  (5.I)  is 


H  =  H^  +  H. ,       IL.  =  eZ  q  +a   , 
-   -0   -1        -c     a^a  ^  ' 


H.  =  m  -^L 


—1 


hz  ^..J    „  ,6  „   .,a  "'"a  + 


(5.5) 


■PHX  apM\  a+p,M-+\^a  %   %% 


where  the  collective  operators  q^^  are  defined  by 

'^a  =  M"2z;^exp(i2rtan/M)qr  1  ,  (5A) 

as  in  Sec.  2.1,  and  obey 

I:V^p];^=0-•    K'%^.=  %-  (5.5) 

For  each  of  our  model  types  (ladder,  ring,  random- coupling,  or 
Hartree-Fock)  the  <I'^g^  are  precisely  the  same  parameters  as  in 
Sec.  2.2,   It  may  be  noted  that  H.  in  (5.5)  bears  a  close  formal 
resemblance  to  (2.l8). 

We  define  the  Heisenberg  operators  q  by 

^Jt)  =  exp[i(t-tQ)H]q^  exp[-i(t-tQ)H)], 

q^"^(t)  =  exp[i(t-tQ)H]qJexp[-i(t-tQ)H]  . 

The  Heisenberg  equation  of  motion  for  any  operator  B  is 
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dB/dt  =  -i[B,  h]_  .  (5-7) 

For  either  the  fermion  or  boson  case  this  yields 
(d/dt  +  iG)q^(t)  =   -iM"^/Ljt)  , 


(5.8) 


where  we  have  used  (5-5)  and  (2.19). 

In  addition  to  the  q  themselves,  we  introduce  the  retarded 
response  or  Green's  operators  G  (t,t'),  defined  by 


Sqjt)  =  Jdt'ZyG^^(t,f)5f^(t')  , 


(5.9) 


G  (t,t')  =0      (t  <  t') 


Here  5f  is  an  arbitrary  infinitesimal  forcing  operator  added  to 

-iM'  \/L  for  times  >  t^,  and  5ci  is  the  increment  in  q  produced 
r  0        a  a 

by  the  addition.   We  restrict  5f  (t)  to  operators  which  anticommute 

(commute)  with  all  the  q  and  q   in  the  fermion  (boson)  ca.se. 
a      a 

The  5f  may  be  regarded  as  arising  from  an  infinitesimal 

7 
modification  of  H.   In  this  way  one  obtains  the  relation' 

G^y(t,f  )  =  [q^(t),  q^^f  )]_^      (t  S  t')  .        (5.10) 

From  (5.6),  (5.10),  and  the  stated  commutation  property  of  8f  , 

it  follows  that  5f  and  6f   commute  with  G   for  either  the  fermion 

7  7  aP 

or  boson  case,  whatever  the  indices  or  time- arguments.   Using  this 
property,  we  obtain 
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(a/at  +  ie)G^^(t,t')   =   -±u\M^^{t,t')  (t  g   f  )    , 

(5.11) 
ir^*^^')  -  Vx%Mx\+3,u-fX  [V'(t,f)$^(t)a,(t) 

+q3"^(t)G^/t,f  )q^(t)  +  qp"^(t)^(t)G^^(t,f  )]  , 


We  shall  also  need  (5-8)  and  (5.11)  in  the  integral  forms 
t 

i(t)  =  5j°^(t)  -  ±M-\  rdfG^°)(t,t')£  (f )         (5.15) 


G^^(t,f)   =   5^^G^°)(t,f)    -   ±M-\j\t"  G^°^(t,t")M^^(t",f).         (5.1^) 

t' 

Here   q  ^    '^(t)    and  5     G^    "^(t^f)    are  the   Heisenberg  destruction  and 

Green's   operators   for  the    'unperturbed'    case  V  =  0.        They  are   given  by 

qj°^(t)   =.  exp[-ie(t-tQ)]q^  (5-15) 

and 

G^°^(t,t')   =     exp[-i6(t-t')]  (t  g   t')    , 

(5.16) 
=0  (t  <  t')    . 


5.5     Correlation  and  Green's   Functions   for  the  Oscillator  Models 
Let  us  now  consider  ensemble   averages   of  the   form 

<fi>H   Tr[ra}    /Tr    (w)    ,  (5.1?) 


^ 


where  B  is  any  operator  and  the  weighting  operator  W  is  a  ftmction 
of  the  Schrodinger  operators  q  =  q.  (*„)•   A  choice  of  W  represents 
a  choice  of  initial  statistical  ensemble  for  the  collection  of 
oscillators.   For  our  present  purposes,  we  restrict  W  to  the 


Gaussian  form 


W  =  exp(-aE  q  ''"q  )  ,  (5-l8) 


where  a  is  a  real  constant  which,  for  bosons  only,  must  further 
obey  a  >  0.   Since  W  then  commutes  with  H,  such  an  initial  ensemble 
is  actually  an  eq\iilibrium  ensemble.   However,  we  shall  not  use 
this  fact  in  developing  the  time -dependent  diagram  analysis.   The 
analogous  Gaussian  weighting  operators  for  the  many-body  problem, 
which  we  shall  introduce  in  Sec.  '^.h,    do  not  commute  with  the  Hamil- 
tonian;  they  represent  genuine  non-equilibrium  ensembles. 

We  define  the  correlation  functions  Q  (t,t')  and  Green's 
functions  G  (t,t')  by 

Q^(t,f)  =  <q;^"^(t')q^(t))  (5.19) 


and 


}Jt,f)  =  <^^(t,f))   .  (5.20) 


We  shall  also  use  the  auxiliary  correlation  functions  Q  (t,t') 
defined  by 

Q^+(t,f)  =  <;^^(t)q^+(f)).  (5.21) 

By  (5-10),  we  have 

G^(t,t')  =  Q^'^(t,f)  +  Qjt,f)     (t  g  f)  .    (5.22) 


55 


It  is  easily  seen  that 

Q  (t,f)  =  Q*{t',t),  Q^^(t,f)  =  Q  +*(t',t)  ,     (5.25) 


and 


C^^^^')  =  ^J^^'^^'P'  (5.2^) 


where  *  denotes  complex  conjugate.   We  see  that  Q  "*"  is  wholly 

determined  by  Q  and  G  .   It  should  be  noted  that  the  time-ordering 
a      a  ° 

operator  T  has  not  been  used  in  any  of  the  definitions  above.   We 
have  not  found  time-ordering  to  be  very  useful  in  the  non-eqioilibrium 
formalism. 

For  V  =  0,  the  correlation  and  Green's  functions  are  independent 
of  a  for  any  M.   They  are  given  by  Q^'^htjV)    and  G^^^t^t'),  where 

Q^°^(t,f)  =  (e%l)-\xp[-i€(t-t')],  (5.25) 

Q+^°)(t,t')  =  [1  +(e%l)"^J  exp[-i€(t-f)]  (5.26) 

and  G^°^(t,t')  is  defined  by  (5.16). 

For  V  ^  0^  G  (t,t')  and  Q  (t,t')  satisfy  the  equations  of  motion, 
and  initial  conditions, 

(S/at  +  ie)G  (t,t')  =  J  (t,t')  ,  (5-27) 

a  a 

G  (tSf)  =  1  ,  (5-28) 

a 

(S/St  +  ie)Q  (t,f)  =  K  (t,t')  ,  (5.29) 

a  a 

Q  (t^^tj  =  (e%l)"^  ,  (5.50) 


where 

+qp"^(t)G^Jt,f  )q^(t)  +  q/(t)q^(t)G^Jt,f  )])        (5.31) 
and 
Kjt,t')   =  -lM-X,x\.3,.+X*a3.X  (C(*')5p(t)a^(t)3,(t)>  .   (5-52) 

These  equations  follow  directly  from  (5'8)  and  (5. 11).   Because  of 
(5.23),  the  equation  for  hGl^{t,t' )/'dt'    is  redundant  with  (5-29). 
The  irreducible  diagram  expansions  for  Jr,('t,t')  and  K^(t,t')  tiim 
out  to  have  a  simpler  form  than  those  for  G(^(t,t')  and  QQ^(t,t')  them- 
selves, and  it  is  for  this  reason  that  we  introduce  the  differential 
equations  (5.27)  and  (5-29)- 

5.4.   Correlation  and  Green's  Functions  for  the  Many-Body  Models 

Let  us  now  return  to  the  many-body  problem.   The  Heisenberg 
equation  of  motion  corresponding  to  the  Hamiltonian  (2.l6),  (2.l8)  is 


t.   -  (5-53) 


(d/dt  +  ±e^)^^^  =  -iM'  ^p^^p^3  -k-s-k+p^r+s 

^   \pnxS3  V^sX  ' 

where  we  have  used  {2.k)    and  (2.19).   In  analogy  to  the  nonlinear 

oscillator  case,  we  define  the  Green's  operators  ^.t^qC"'^''''' )  ^y 
t 
Sq^Jt)  =/dfE^p\^^p(t,t')6f^p(t'),  (5.34) 
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where  Sf^^^  is  an  infinitesimal  forcing  operator  introduced  on  the 
right  side  of  (5.35). 

We  again  consider  ensemble  averages  of  the  form  (5.1?),   where 
now  the  weighting  operator  has  the  form 

^=  ^^Ka^(^^\J  '  (5.35) 

with 

We  restrict  the   real  function  w  (k)  by  the   condition 

w(k)   =    |w(k)|   cc      e^  (k->  »)    ,  (5.56) 

For  bosons   only,    we   further  require  w(k)    >  0  for  all  k. 

The   correlation  and  Green's   functions   for  the  many-body  models  may 

now  be  defined  by 

W^'^'^  =<^ka'^*')\J*>^       \>'^')=   <\Jt)5j(f)>,     (5.37) 

\J^'*')   =   <fkkaa(^^*'i>    •  (5.58) 

For  V(x)  =  0,  the  correlation  and  Green's  functions  are  inde- 
pendent of  a  for  any  M;,  as  was  the  case  for  the  nonlinear  oscillator. 
They  are  given  by 

Sc^°^(t,t' )  =  N^(°^exp[-ie^(t-t-)],  (5.59) 

Qj^+(°^(t,t')  =  (1  +  N^(0))  exp[-iej^(t-t')]  (5.1,0) 
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jj°^(t,t')  =  exp[-iejjt-f)]       (t  g  f  )  , 
=0  (t  <  t')  , 


where 


(0) 


=   (e^^^V  1) 


For  V(x)  ^  0,  we  have 


(B/St  +  iei,)\^(t,t' )  =  Jj,^(t.t') 


G^JtSf)  =  1, 


(5^1) 


(5^2) 


(5.^5) 


(S/St  +  ie^)C^_^(t,t' )  =  Kj^^(t,f  ). 


WV^o)  =  \ 


(0) 


(5.H) 


where  J  (t,t' )  and  K,  (t,t')  are  given  by  obvious  modifications  of 
(5.31)  and  (5.32).  The  analogs  of  (5.22)-(5.2i|-)  hold  also  for  the 
many-body  case. 

The  partition  of  particles  and  of  energy  among  the  momentum  modes 
as  a  fimction  of  time  may  be  expressed  directly  in  terms  of  ft  (t,t'). 
We  have,  immediately,, 


{l)=W\J^'^^ 


(5.^5) 


In  analogy  to  IiC+.lO),  we  may  easily  establish  from  (5-33)  the 
relation 

ao^jt.t-)- 


{^) -K\  K\J^^^^  ^ 


^F 


(5.^6) 
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5.6  Approach  to  Equilibrium 

Let  us  consider  the  case  t  -»•  -»,   Then  at  finite  t  ve 
anticipate  that  the  ensemble  specified  by  (5-55)  will  have  evolved 
into  a  state  of  statistical  stationarity.    This  need  not  be  so  for 
the  random- coupling  model,  which  may  be  unstable  to  catastrophic 
collapse  because  of  the  unboundedness  of  H. .   However,  it  seems 
assured  for  the  ladder  and  ring  models,  if  V(x)  satisfies  the  con- 
ditions which  make  H.  bounded  from  below.  If  an  approach  to 
equilibrium  is  granted,  we  anticipate  for  any  finite  t  and  t'  that 
^^(t,t')  and  \a_(t,,t')  will  take  the  forms 

It  should  be  noted  that  we  are  not  invoking  an  adiabatic  switch-on 
of  V(x).   We  are  simply  fixing  the  initial  statistical  ensemble  by 
a  choice  of  W  and  then  letting  the  ensemble  evolve  according  to  the 
exact  equations  of  motion. 

The  eqmlibriuni  ensemble  which  is  achieved  will  in  general 
not  be  a  grand  canonical  ensemble.    The  reason  is  that  (5.55)  does 
not  in  general  represent  an  initial  distribution  of  total  energy  (kinetic 
plus  potential)  and  total  particle-number  corresponding  to  a  grand 
canonical  ensemble.   Since  H  and  N  are  constants  of  motion,  the 
distribution  of  these  quantities  can  not  change  with  time.   Never- 
theless, we  conjecture  that  for  our  models  in  the  limit  M  -►  »  the 
equilibrium  functions  Q,  (t-f)  and  G^q,^*"^'^  ^"^"'"  ^^  identical  in 
value  with  corresponding  averages  defined  over  the  grand  canonical 
ensemble  with  the  same  /H")  and  /N^  .   In  the  limit  M  ->■  »,  both  the 
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grand  canonical  ensemble  and  the  ensemble  specified  "by  (5*55) 
represent  total'-energy  and  total-particle-number  distributions 
which  are  peaked  with  infinite  sharpness  about  their  means.   We  have 
noted  previously  [cf.  (2.27)  and  (2.28)]  that  the  M  systems  in  the 
collection  are  coupled  for  our  models;  they  exchange  both  energy 
and  particles.    This  means  that,  with  respect  to  the  achievement 
of  statistical  equilibrium,  the  models  represent  infinitely  large 
super-systems  in  the  limit  M  -*  «>,  even  if  the  mean  nimiber  of  particles 
per  system  is  small.   It  is  on  the  basis  of  these  facts  that  we 
conjecture  the  identity  of  %^'^^-'^')    and  G^  (t-f)  for  the  two  kinds 
of  ensemble.   The  conject\i3re  supposes  certain  ergodic  properties, 
as  does  the  assumption  that  ensembles  specified  by  (5-35)  will  reach 
equilibrium  at  all. 

If  we  turn  from  the  models  to  the  true  problem  (all  '*^^j3lJ^=  1)> 
then  the  M  systems  in  the  collection  are  not  coupled,  and,  if  the 
mean  number  of  particles  per  system  is  small,  we  do  not  have  a 
large  super-system  in  any  dynamical  sense.   In  this  case,  the  grand 
canonical  ensemble  and  the  ensemble  specified  by  (5-55)  cannot  be 
expected  to  yield  the  same  averages  when  t  -»  -<».   It  does  not  follow 
that  our  Gaussian  ensembles  are  physically  inappropriate.   The  grand 
canonical  ensemble  is  used  for  equilibrium  calculations  more  because 
it  is  mathematically  convenient  than  because  it  is  -uniquely  appro- 
priate physically.   For  our  present  purposes,  Gaussian  ensembles 
are  the  ones  which  are  most  convenient  mathematically.   The  real 
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justification  for  either  choice  of  ensemble  is  the  hope  that  for  a 
dynamically  large  system  (M~  /n\-»  oo  for  the  true  problem;  M  -*  »  and 
M  (^y   finite  for  the  models)  the  values  of  physically  interesting 
equilibrium  averages  are  insensitive  to  a  substantial  range  of  choices 
of  ensemble.   We  shall  not  attempt  to  go  further  into  this  matter 
here.    The  questions  which  arise  are  not  unique  to  our  investigation. 

The  conjecture  made  above  concerning  O.  (t-f)  and  G,  (t-f) 
for  the  models  implies  that  these  quantities  may  be  identified  with 
equilibrium  correlation  and  Green's  functions  of  types  discussed 
previously  by  a  number  of  authors.  The  temperature -domain  propa- 
gators S,  (u,u')  may  be  expressed  in  terms  of  the  equilibrium  cor- 
relation and  Green's  functions  by  analytic  continuation  of  the  latter. 
The  continuation  is  given  formally  by  the  relations 

^°^  ^"^  (5.48) 

=  +\J-i^  +  iu')        (u  S  u')  , 

which  follow  from  (j-'^)  and  (5.57)  if  (^   \   is  taken  throughout  as  an 
average  over  the  grand  canonical  ensemble. 

There  is  a  useful  general  relation  between  the  equilibrium 
functions  0_  (t-f)  and  G,  (t-f).   If  we  define  G.  ^(t-f )  by 

=  Gj^J(t'-t)      (t  <  f  )  , 


k2 


then   the   transforms 


K^^""^      =      (2rt)-yG^^^(t)exp(lcot)dt    , 


(5.50) 


,(-)        =      (2«)  Q^Jt)exp(i.. 


it)dt 


are   related,  by 


(.)        =      (eP^-^).   l)"\a'^-)    •  ^5.51) 


Equation  (5.5I)  has  just  the  form  of  the  free-particle  Fermi-Dirac 

~  S 
or  Einstein-Bose  distribution  law  provided  that  we  interpret  G   (w) 

as  a  density  of  states  for  momentum  k  and  '^~(^)  as  a  mean  occupancy. 

For  the  grand  canonical  ensemble,  (5. 51)  may  be  obtained  directly 

from  the  definitions  of  G,   (w)  and  Qj^q^('^')  "by  using  the  cyclic 

11  12 
properties  of  the  trace.  '  However,  (5.51)  may  also  be  deduced 

as  a  necessary  condition  for  equilibrium  under  coupling  to  a  ther- 
mometer, without  specifying  the  precise  nature  of  the  equilibrium 
ensemble.    This  is  clearly  preferable  for  the  present  application. 
We  hope  to  present  a  derivation  of  this  type  at  another  time.   The 
nature  of  the  argument  has  been  stated  previously. 


6.   DIAGRAM  EXPANSIONS  FOR  THE  NONEQUILIBRIUM  CORRELATION  AND 

GREEN'S  FUNCTIONS 
6.1.   Primitive  Diagram  Expansions 

The  generation  of  primitive  and  irreducible  linked-diagram 
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expansions  for  ^-^Ji^^^')    and  K:^^(t,t')  is  straightforward  in 
principle,  but  somewhat  intricate  in  practice.   We  shall  introduce 
•the  procedure  by  means  of  the  nonlinear  oscillator  example  dis- 
cussed in  Sec.  5.2.   Suppose  that  we  carry  out  a  formal  iteration 
solution  of  the  integral  equations  (5. 15)  and  (5.l4).   We  thereby 
obtain  q  (t)  and  G  (t^f)  expressed  as  power-series  in  V.   The 
coefficients  in  these  series  consist  of  terms  of  the  following 
kind:   Each  term  is  a  multiple  integral  over  a  product  of  Kr(Snecker 
s;yTnbols,  of  factors  of  the  form  M  ^,    and  of  unperturbed  operators 


-  ^  ,     .    -        ,  --  -    -   The  term  is  summed 


of  the  form  q  S^\    ^  /^^^  0^°),  or  0^°)^ 
over  the  intermediate  indices.  Jf  the  power-series  are  substituted 
into  (5- 51)  and  (5.52),  we  obtain  expressions  for  J  (t,t')  and 
K  (t,t')  as  power-series  in  V.   The  terms  which  make  up  the  coef- 
ficients in  these  series  involve  integrals  over  ensemble-averages 
of  products  of  factors  q_'  ,^^\    q.  ,^^^\    G^'^\    and  G^^^*. 

In  order  to  evaluate  the  ensemble-averages^  we  first  note 
that  G^  '^(±,t'),  given  by  (5.I6),  is  a  c-number  function  and  may 
be  taken  outside  the  brackets  \  /  .   The  remaining  averages  over 
q^   and  q     factors  may  be  evaluated  by  using  a  statistical  form 
of  Wick's  theorem,  which  follows  from  the  Gaussian  form  of  the 
weighting  operator  (5.I8).    We  pair  the  q^    and  q"*"^  '  factors 
in  all  possible  ways^  maintaining  always  the  original  left-right 
order  of  the  two  factors  in  a  pair.   We  replace  each  pair  by  its 
individual  ensemble-average  of  the  form 
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<a,^(°)(t.)a/°)(t)>  =  6^Q^°)(t,t.) .  (6.1) 

<;a^(o)(t)a,^^°)(f))  =  6^^Q-^^°)(t,f) ,  (6.2) 

or 

<q;°'(t)qV°)(f))    =   <S;'°'{t)4/(«){t.)>     =   0  .      (6.3) 

Then  we  take  the  product  of  all  the  individual  averages  for  each 

pairing  and  sum  over  all  the  possible  pairings.   For  fermions  only, 

P 
we  multiply  each  product,  before  summing,  by  (-1)  ,  where  P  is  the 

number  of  permutations  required  to  obtain  the  particular  pairing. 

Let  us  suppose  that  we  have  carried  out  the  procedure  just 
described  and  have  then  performed  the  summation  over  all  the  inter- 
mediate indices  which  arise.   Each  contribution  proportional  to  V  , 
in  the  expansion  for  either  J  (t,t')  or  K^(t,t'),  then  consists  of 
a  multiple  integral  over  a  product  of  Q^   ,  Q  ^   ,  G^   ,  and  G^ 
functions,  multiplied  by  a  stmimed  product  of  KrSnecker  symbols  and 
n  factors  m"  $.    The  Kronecker  symbols  express  the  'collective -index 
conservation'  which  also  characterized  the  equilibrium  analysis  of 
Sec.  5.   A  consequence  is  that  the  sums  over  products  of  Rrftliecker 
symbols  and  factors  m"  $  turn  out  to  be  precisely  the  quantities 
C  ,  (a)  of  the  equilibriim  theory.    These  quantities  contain  the 
entire  dependence  of  the  contributions  upon  the  $'s  and  upon  a. 

As  the  preceding  paragraph  suggests,  the  formal  expansions 
for  J  (t,t')  and  K  (t,t')  in  powers  of  V  have  systematic  diagram 
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representations  which  resemble  the  primitive  linked- diagram  expan- 
sion of  the  equilibriijm  theory.   It  is  immediately  apparent  that 
only  linked-diagraras  arise.   Any  intermediate  index  which  occurs 
arises  from  an  iterative  branching  of  the  Heisenberg  equations  of 
motion  and  thus  is  necessarily  linked  by  a  Kronecker  symbol  to 
indices  which  occurred  previously  in  the  iteration  process. 

In  order  to  write  the  complete  contribution  to  J  (t^f)  or 
K  (t,t')  which  is  proportional  to  V^,  we  first  write  down  all  the 
distinct  nth-order  primitive  linked  diagrams,  just  as  in  the  equil- 
ibri\jm  treatment.   Now,  however,  it  turns  out  that  each  diagram 
Hand  hence  each  C   (a)1  above  the  first  order  is  associated  with 
more  than  one  contribution.   The  method  of  forming  the  diagram 
contributions  is  most  clearly  indicated  by  giving  some  examples. 
Consider  the  first-order  contribution  to  K  (t,t')  associated  with 
the  diagram  of  Fig.  ^.   The  one  vertex  in  the  diagram  is  associated 
with  the  $  factor  that  appears  explicitly  in  (5.32).   The  contribution 
is  obtained  from  the  product  of  the  zeroth-order  terms  in  the 
iteration  expansions  of  all  the  operators  in  (5.52).   Tliis  is 
because  a  factor  V  already  appears  explicitly  in  (5-32).   There  are 
two  possible  non-vanishing  pairings  of  the  creation  and  destruction 
operators  according  to  Wick's  theorem: 

(aX)(pji)    and    (a|i)(3X)  . 

Only  the  first  pairing  corresponds  to  the  present  diagram.   By  (6.1), 
it  yields  the  contribution 
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-iC^.^(a)VQ^°)(t,t)Q^°)(t.t')  ,  {G.h) 

where  C    (a)  is  given  by  (3-15) •   The  second  pairing  corresponds  to 
the  exchange  diagram  Fig.  5?  and  involves  the  quantity 

We  have  diagrammed  the  contribution  (6.'+)  by  labeling  the 
vertex  in  Fig-  ■^  with  time  t,    labeling  the  outer  end  of  the  incoming 
external  line  with  time  t',  and  writing  the  two  Qc        factors  along 
the  lines  with  which  they  are  associated.   No  factor  is  associated 
with  the  outgoing  external  line,  and  tliis  is  true  of  every  contri- 
bution to  the  expansions  of  K  (t,t')  and  J  (t,t').   The  reason  is 
that  we  are  dealing  with  the  differential  equations  (5.29)  and  (5.27) 
instead  of  with  Q  (t,t')  and  G  (t,t')  directly.   In  the  box  on  the 
right  side  of  Fig.  ^,  we  have  given  an  alternative  s;>Tnbolic  represen- 
tation of  the  contribution  of  this  diagram  to  (5.29). 

In  Fig.  6,  we  have  diagrammed  the  contribution  to  J  (t,t') 
associated  with  C    (a)  •   Tliis  contribution  arises  from  the  third 
term  on  the  right  side  of  (5-3l)>  and  it  is  obtained  by  replacing  all 
the  operators  in  that  term  by  their  zeroth-order  values.    By  (6.1);, 
the  contribution  therefore  is 

-lC^.^(a)VQ^°^(t,t)G^°)(t,f)  .  (6.6) 

A  symbolic  representation  of  the  contribution  of  this  diagram  to 
(5.27)  is  given  in  the  box  on  the  right  side  of  Fig.  6.   (The  under- 
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lined  a  represents  G  .)   The  second  term  on  the  right  side  of  (5-51) 
gives  rise  to  the  first-order  exchange  contribution  to  J  (t^f),  associated 
with  Fig.  5  and  C,  „(a).   By  (6.5);.  the  first  term  on  the  right  side 
of  (5«5l)  gives  no  first-order  contribution. 

There  are  three  contributions  to  J  (t.f)  associated  with  the 
second-order  diagram  shown  in  Fig.  7a-   They  are  diagrammed  in  Figs. 
7b,  7c,  and  7(1  •   In  the  boxes  beside  the  diagrams  we  have  represented 
symbolically  the  particular  iteration-substitutions  (iterative  branch- 
ings) and  operator-pairings  associated  with  the  contributions.   Below 
this,  we  have  shown  (in  square  brackets)  the  left-right  ordering  of 
the  final  set  of  creation  and  destruction  operators  which  are  pro- 
duced, in  each  case,  by  the  iteration-substitutions.   This  must  be 
kept  account  of  to  determine  whether  Q^  "^  or  Q  ^    factors  are  associ- 
ated with  given  lines  and  to  determine  the  sign  of  the  contribution 
in  the  fermion  case.   The  total  contribution  of  Fig.  7  to  J  (t,t')  is 

C2.,(a)v2|  [-q(0)(,.,  ,t)Q+(0)(t,t")G^°)(t.t-  ) 
t' 

lQ(°)(t",t)G(°)(t,t")Q^°^(t.t") 

:iiG^°)*(t,t")  Q^°^(t,t")  Q^°\t,t")]G(°^(t",t')'if,    (6.7) 
where 

The  three  terms  which  comprise  (6.7)  are  associated,  from  left  to 
right,  with  Figs.  7b,  7c,  and  7d.,  respectively. 


The  contributions  to  K  (t,t')  and  J  (t,t')  associated  with 
all  the  primitive  diagrams  may  be  determined  in  the  fashion  illus- 
trated above.   For  each  diagram^  one  traces  through  all  the  iter- 
ative branchings  and  non-vanishing  operator-pairings  which  correspond 
to  the  diagram  topology.   Then  one  writes  down  the  contributions  by 
using  (6.1)  and  (6.2).   The  primitive-diagram  contributions  have  the 
following  general  characteristics.   Each  nth  order  contribution  to 

J  (t,t')  contains  a  factor  C  _  (a)  which  multiplies  an  n-fold  integral 
a  n;p 

over  n  factors  G    or  G     (in  some  combination)  and  n  factors  Q 

or  Q    .   The  contributions  to  K  (t,t')  are  similar  except  that  there 

are  n-1  factors  G    or  G     and  n+1  factors  Q    or  Q    .   In  every 

contribution  to  J  (t,t')^  but  in  only  some  of  the  contributions  to 

K  (t,t')^  the  factor  associated  with  the  incoming  external  line  is  a 

G    factor.   We  have  already  remarked  that  no  factor  is  ever  associated 

with  the  outgoing  external  line  in  either  the  J  (t.t')  or  K  (t.t') 

a  a 

expansions. 

The  primitive  linked-diagram  expansions  for  the  fimctions 
J  (t,t')  and  K  (t,t'),  which  appear  in  (5.i^-3)  and  (5-^^)^  may  be 
obtained  in  close  analogy  to  the  analysis  for  the  nonlinear  oscillator. 
The  only  difference  is  that  there  is  a  momentum  associated  with  each 
line  in  the  diagrams,  and  the  results  must  be  summed  over  the  inter- 
mediate momenta,  as  in  the  equilibrium  analysis.   For  example,  the 
contributions  to  J  (t,t')  which  correspond  to  (6.6)    and  to  the  first 
term  in  (6.7)  are 

-iC^.^(a)ZpVQQ^^°)  it,t)Gj^^  (t,t  •)  (6.9) 
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X  Gg^°\t,f  >  )G^^°\f  S  f  )df  •  ,      (6.10) 

respectively. 

Let  us  now  consider  the  limit  M  -♦  »  .   As  we  discussed  in 

Sec.  k,    all  the  primitive  linked  diagrams  may  be  divided  into  two 

classes,  'wanted'  and  'unwanted, '  for  each  of  the  stochastic  models. 

The  C  .  (a)  for  wanted  diagrams  all  have  the  value  one.   For  any 

unwanted  diagram  of  finite  order  n,  we  have  C   (a)  -»■  0  as  M  -♦  «  . 

n;p 

Let  us  extend  to  the  non-equilibri\im  case  the  fundamental  assumption 
made  in  Sec.  5  about  unwanted  diagrams.  We  assume  that  the  total 
contribution  to  J  (t,t'),  K^(t,t'),  Jj^ctC't;.  t '  ) ,  and  \J.^,^')    from 
all  unwanted  diagrams  vanishes  in  the  limit  M  -»  «  .  We  shall  reserve 
all  discussion  of  the  validity  of  this  assumption  for  Sec.  7.  An 
immediate  implication  of  the  assumption  is  that  J  (t,t'),  K  (t,t'), 
G  (t^f),  Q  (t,t'),  and  the  corresponding  functions  for  the  many- 
body  problem  all  become  independent  of  a  in  the  limit.   We  may 
therefore  omit  the  index  a  in  these  functions.   Then,  in  analogy 
to  (5.21),  we  find 


<«[n]'(^')^[m](^)> 


(6.11) 


(^[m](*)^[n]'(^'^   =^rS^^'^')> 


with  corresponding  relations  for  the  many-body  problem. 
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6.2  Irreducible  Diagram  Expansions 

Irreducible  linked-diagram  expansions  for  J^(t,t')  and 
K  (t,t')  may  be  constructed  by  the  following  rules:   Retain  only 
the  irreducible  diagrams;  that  is,  those  without  self-energy  parts. 
(See  Rule  2'  of  I:  Sec.  h   for  the  definition  of  an  irreducible 
diagram  and  of  a  self-energy  part.)   Then  alter  the  primitive 
contributions  associated  with  these  diagrams  by  replacing  each  factor 
q(0)^  (.(0)*^  ^(0)^  ^^  ^+(0)  ^j^g^ein  ^^^   a.  factor  G^,  G^,  Q^,  or  q/ 
having  the  same  time -arguments.   Here  a   is  the  collective  index  that 
labels  the  line  associated  with  the  factor.   [The  external  line 
(a  =  a)  does  not  form  an  exception  to  this  rule  .J    Corresponding 
rules  hold  for  the  irreducible  expansions  for  J   (t,t')  and 
Vt,t'). 

The  formal  validity  of  these  irreducible-diagram  expansions 

may  be  demonstrated  in  two  ways,  neither  of  which  will  be  carried 

out  here.    The  first  way  is  to  construct  directly,  by  iteration 

and  the  use  of  Wick's  theorem,  the  primitive-diagram  expansions  for 

G  (t.f),  Q  (t,t'),  etc.     One  then  substitutes  these  expansions 
0  a 

into  the  irreducible  expansions  for  J  (t,t')  and  K  (t,t')  and  compares 

the  results  with  the  primitive  expansions  for  J^(t,t')  and  K^(t,t'). 

The  second  method  does  not  require  the  explicit  primitive 

expansions  for  G  (t,t')  and  Q  (t,t').   It  involves  calculating 

the  changes  in  J  (t,t')  and  K  (t,t')  produced  by  small  variations 

in  the  $    ,  and  it  gives  directly  the  contribution  proportional  to 
apilX 
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each  irreducible  C^   (a).   The  method  is  illustrated,  for  another 

17 
application,  in  reference  2. 

The  irreducible  expansions  are  formally  exact  for  any  M.   In 
the  limit  M  -*  «=,  our  assumption  that  unwanted  diagrams  make  no  con- 
tribution to  the  primitive  expansion  implies  that  these  diagrams  make 
no  contribution  to  the  irreducible  expansions  also.   The  argiment  is 
precisely  analogous  to  that  given  in  Sec.  5  for  the  equilibrium  case. 
Thus,  we  obtain  the  formal  irreducible  expansions  for  J  (t,t')  and 
K  (t,t')  in  the  limit  M  ^  »  by  retaining  only  the  wanted  irreducible 
diagrams.   For  each  of  the  stochastic  models,  these  expansions  can  be 
summed  to  give  closed  equations  for  G(t,t')  and  Q(t,t'),  in  analogy 
to  the  equilibrium  case.   The  results  are  substantially  more  compli- 
cated than  in  the  equilibriimi  case,  however.   The  complete  equations 
for  all  the  models  are  given  in  Appendix  B,  for  both  the  nonlinear 
oscillator  example  and  the  many -body  problem. 

We  have  remarked  earlier  that  the  irreducible  expansions  for 

J  (t,t')  and  K  (t,t')  are  simpler  than  those  for  G  (t,t')  and  Q  (t,t') 

themselves.   The  expansion  for  G  (t,t')  actually  may  be  constructed 

very  easily  from  that  for  J  (t,t')  by  integrating  (5.27).   Thus  we  find 

t 
Gjt,f)  =  G^°^(t,t' )  +   /   G^°^(t,t")  Jjt"  ,t')dt"  .    (6.12) 
t' 

The  factor  G^  '^(t,t")  in  the  integrand  may  be  considered  to  be  associ- 
ated with  the  outgoing  external  line  in  each  of  the  irreducible  dia- 
grams.  [We  recall  that  no  factor  was  associated  with  this  line  in  any 
of  the  diagram-contributions  to  J  (t,t').]   However,  we  have  not  sue- 
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ceeded  in  finding  an  equally  compact  irreducible  expansion  for  Q  (t^f), 
Equation  (5.29)  gives  a  substantially  more  complicated  result  upon 
integration  than  does  (5.27). 

The  difficulty  in  constructing  a  compact  irreducible  expan- 
sion for  Q  (t,t')  is  apparent  from  the  structure  of  the  primitive 
expansion  for  this  quantity.   In  the  primitive  expansion  for  G  (t^f), 
the  factors  associated  with  the  incoming  and  outgoing  external  lines 
are  always  G^   factors.   This  property  underlies  (6.12)  [which  is 
analogous  to  I:  (4.12)].   In  the  primitive  ejcpansion  for  Q  (t^f), 
some  contributions  have  a  G    factor  associated  with  the  outgoing 

fo)-^ 

external  line  and  a  G^     factor  associated  with  the  incoming  external 
line.   Tlie  remaining  contributions  have  a  Q^   factor  associated 
with  either  the  incoming  or  outgoing  external  line  and  a  G 
or  G    factor  associated  with  the  other  external  line.   This  precludes 
an  expression  for  Q  (t,t')  of  the  simple  form  (6.12). 

It  is  clear  from  what  has  been  presented  in  this  Section  that 
our  non-equilibrium  primitive  and  irreducible  expansions  do  not  have 
the  simplicity  and  compactness  of  the  diagram  analysis  for  the  equil- 
ibrium case,  which  we  discussed  in  Sec.  3'    This  may  be  an  unavoidable 
penalty  for  abandoning  equilibrium,  but  it  may  also  be  that  our  formu- 

-|  o 

lation  is  unnecessarily  awkward. 

7.   VALIDITY  OF  THE  CLOSED  MODEL  Et^UATIONS 

7.1  Description  of  Method 

We  wish  in  Sec.  7  to  investigate  the  basic  assumption  that  the 
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total  contribution  of  the  ;inwanted  diagrams  vanishes  in  the  limit 
M  -*  00  .   This  assumption  was  stated  for  the  equilibrium  case  in 
Sec.  5  and  for  non-equilibriiun  in  Sec.  6.   It  was  the  essential 
ingredient  in  establishing  closed  equations  for  the  model  propa- 
gators, correlation  functions,  and  Green's  functions.   The  procedure 
we  shall  use  here  is  to  integrate  the  Heisenberg  equations  of 
motion  by  replacing  them  with  a  set  of  difference  equations  in- 
volving the  discrete  times  t  ,  tQ  +  At,  t^  +  a^t,  ....  At  each 
stage  of  integration,  this  finite-difference  method  yields  approxima- 
tions to  the  non-equilibrium  correlation  and  Green's  functions 
which  contain  as  coefficients  the  same  quantities  C  .  (cc)  that 

arose  in  the  iteration  procedure  of  Sec.  6.   Only  C  ,  (a)  of 

n;p 

finite  order  appear  after  a  finite  number  of  integration  steps. 

The  finite-difference  method  of  integration  actually  consti- 
tutes a  definition  of  the  Heisenberg  equations  of  motion.   Con- 
sequently, we  may  hope  that  it  converges  in  the  limit  At  ->  0  to 
yield  the  exact  correlation  and  Green's  functions,  whether  or  not 
the  iteration  solutions  of  Sec.  6  converge.   If,  moreover,  the 
convergence  as  At  -►  0  is  independent  of  M  for  large  M,  we  may  con- 
clude immediately  that  the  unwanted  diagrams  make  zero  total  contri- 
bution to  the  exact  correlation  and  Green's  functions  in  the  limit 
M  ->  00.    This  is  because  the  unwanted  C  ,  (a)  which  appear  in  any 
given  finite-difference  approximation  are  of  finite  order  and  vanish 
in  the  limit  M  ^  oo. 

In  what  follows,  we  do  not  attempt  to  prove  rigorously  the 
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requisite  convergence  properties  of  the  finite -difference  approxi- 
mations.  We  make  the  convergence  plausible  by  showing  that  the 
time-derivatives  of  the  matrix-elements  of  the  q   (t)  have  bounds 
which  are  independent  of  M  in  an  appropriate  mean-square  sense. 
This  is  done  for  all  t  without  any  appeal  to  perturbation  expansions. 
However,  it  is  essential  to  our  analysis  that  the  Hamlltonian  be 
modified  by  removing  all  momentum  modes  above  an  arbitrarily  high 

but  finite  cut-off  k    .   The  significance  of  the  momentum  cut-off 
max 

is  discussed  for  the  several  stochastic  models.   In  the  case  of 

the  ring  and  ladder  models,  it  is  concluded  that  the  exact  functions 

G,  (t.f)  and  a  (t,t')  for  given  k  should  be  negligibly  dependent 

on  k    if  k    is  high  enough  and  if  V(x)  satisfies  the  conditions, 
max     max      b      t=  \  / 

stated  in  Sec.  2,  which  yield  lower  bounds  to  the  eigenvalues  of 
the  model  Hamiltonians . 

After  establishing,  to  the  extent  described,  that  the  un- 
wanted diagrams  do  not  contribute  to  the  non-equilibrium  correlation 
and  Green's  functions,  we  point  out  that  this  does  not  complete 
the  justification  of  our  formal  closed  equations  for  the  stochastic 
models.   If  the  wanted  diagrams  are  sufficiently  numerous  in  high 
orders  that  they  form  non- converging  series,  then  a  \iniqueness 
question  arises  in  the  summation  of  these  diagrams.   We  treat  this 
question  by  regarding  the  closed  integro-differential  equations  them- 
selves as  limits  of  finite-difference  equations,  rather  than  as  sum- 
mations of  infinite  classes  of  perturbation  terms.   Finally,  we  examine 
the  conditions  under  which  the  non-equilibrium  closed  model  equations 
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imply  those   for  eqiailibritun. 

7-2     Bounds   on  Matrix  Elements   and  Their  Time-Derivatives 
Let  \is   consider  the  positive -definite  quantities 

(^^Fjt)   =     <qj(t)qjtj>       =      Qjt,t)  (r  =   0)    , 

(7.1) 


a' 


where  the  q  (t)  are  the  Heisenberg  operators  for  the  nonlinear 
oscillator  example.  We  have  noted  previously  that  the  weighting 
operator  (5-18)  commutes  with  the  Hamiltonian.   It  follows  that 


Fjt)      =      ^^^^FjtQ)    (r  =  0,  1,  2,  ...)  (7.2) 


^^^F  (t)   =   ( 
for  all  t.   By  (5.30),  we  then  have 


^°^Fjt)  =  N  ,  (7.5) 

where  N  =   (e^  +  1)""^.      Using   (7.2),    (5.6),    (5.8),    and  Wick's 
theorem,   we   find 

^l^Fjt)   =  1  (e   +  VN    [C^.^(a)    -    C^.2(a)j}  ^ 

(7.M 
V%^(1  +  N)[C2.-L(a)    +  C^.^icL)]    . 

The  C   (a)  which  appear  in  (7-^)  have  been  defined  by  (5. 15), 
(3.17),  (6.5),  and  (6.8).   It  follows  from  (2.19)  that  they  are  all 
real. 
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(r) 

Expressions  similar  to  (7.^)  may  be  foimd  for  the  ^  ^F  ft) 

(t) 
of  any  finite  order  r.   Each    F  .ft)  may  be  evaluated  for  all  t 

by  repeated  differentiation  and  self-substitution  of  ('5-8)  at  t  =  t  , 

folloved  by  the  use  of  Wick's  theorem.   Tlie  result  is  a  polynomial  of 

finite  degree  in  1  and  the  C  .  (a).   Only  the  C  ,  (a)  of  order  n  S  2r 

(r) 
appear.   It  follows  from  (4.3)  that  each^  -^F  ft)  has  a  finite  bound 

which  depends  on  N  but  is  independent  of  t,  a,  and  M.,  and  is  also  inde- 
pendent of  the  choice  of  model.   In  particular,  this  is  true  in  the 
limit  M  -♦  00. 

In  order  to  make  clear  the  significance  of  this  result,  let 

fr) 
us  write  the  traces  which  define  the  ^   F  (t)  in  the  explicit  fonn 

(O^Fjt)  =  Z-lE^^^,exp[-aN(s)]  j  ^  '  |qj  t)  |  s>|  ^  , 

^^^Fjt)  =  Z-\^,exp[-aN(s)]  j  ^  ■  |d\(t)/dt^|  s>|  ^  (7.5) 

=  Z-\^,exp[-aW(s)]  ]  d^  <s  '  |  qjt)  |  s  > /dt^|  ^   (r  g  1)  , 

where 

Z  =  Z  exp[-aN(s)]  (7-6) 

and  (s'|q  ft)|s)  is  the  matrix  element  of  q^ft)  between  the  states 
s  and  s'.   The  sums  are  over  the  complete  set  of  joint-eigenstates 
s.s'  of  the  Schrodinger  number-operators  q^  q^,  and 

N(s)=   <s|E^qJciJs>  (7.7) 

is  the  total  number  of  quanta  in  the  state  s. 
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(r) 
The  boimdedness  of  ^   F  (t)  implies  first  of  all  that  the  sums 

in  (7-5)  converge  to  finite  limits.   In  parti cialar,  it  implies  that 
the  Slim  over  the  complete  sets  of  states  s  which  have  successively- 
increasing  values  of  N(s)  converges  as  N(s)  -*  m.      The  density- of- states 
in  the  space  of  the  occupation  numbers  is  such  that  the  factor 
Z  exp[-aN(s)]  in  (7-5)  gives  the  principal  total  weight  in  the  sura 
to  states  s  with  N(s)  =  0(MN)  .   In  the  limit  M  -♦  oo^  the  sum  is  in 
effect  sharply  confined  to  states  s  with  N(s)/m  =  N.   Let  us  consider 

the  complete  set  of  states  with  a  given  value  of  N(s)/m.   The  bounded- 

fr) 
ness  of  ^   F  (t)  then  implies  that  the  rth  time-derivatives  of  the 

complete  set  of  matrix  elements   <^  s'|q  (t)|sN  have  a  mean-square 
bound  which  is  independent  of  M  in  the  limit  M  -♦  <»  and  which  is  finite 
if  N(s)/m  is  finite.   It  should  be  noted  that  the  matrix  elements 
<(_s'|q  (t)|s^  are  non-vanishing  only  if  N(s')  =  N(s)  -  1.   This  fol- 
lows from  (5.6)  and  the  fact  that  N  commutes  with  H. 

We  must  now  extend  the  analysis  to  the  many-body  problem,  for 
which  H  is  given  by  (2.l6)  and  (2.l8).   We  seek  mean-square  bounds  on 
the  time-derivatives  of  the  matrix  elements  '^s'|q  (t)|s_)  ,   where 
s  and  s'  are  now  joint-eigenstates  of  all  the  number-operators 


case  only  if  a  momentiom  cut-off  is  introduced,  as  mentioned  in  Sec.  7-1- 

We  remove  from  H  all  terms  which  contain  any  q  or  q  factor  whose  momen- 

ttmi  index  exceeds  in  magnitude  some  arbitrarily  high  but  finite  value 

k   .   It  is  clear  that  the  cut-off  preserves  the  Hermiticity  of  H. 
max  — 
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We  shall  reserve  for  Sec-  7-5  all  discussion  of  the  dynamical  signifi- 
cance of  the  cut-off  for  the  several  stochastic  models. 

Bounds  for  the  many-body  matrix  elements  with  the  momentum 
cut-off  imposed  can  be  obtained  in  close  analogy  to  the  procedure  for 
the  nonlinear  oscillator.   We  introduce  the  special  Gaussian  weight 
operator 

W-  =  exp[-  a  \^\^]    =  exp[-  oN]  ,  (7.8) 

where  a  is  a  real  constant.   (For  bosons  only,  a   must  also  be  posi- 
tive.)  Unlike  the  more  general  form  (5-55):>  this  operator  commutes 
with  H  and  therefore  corresponds  to  an  equilibrium  ensemble.   It  may 
be  regarded  as  the  infinite -temperature  limit  (p  -♦  0,  -^ii  -*  o)    of  the 
weighting  operator  exp[-p(H  -  hN)]  which  yields  the  grand  canonical 
ensemble  for  the  problem  with  momentum  cut-off. 

Let  (  /"  '  denote  a  trace  weighted  by  W' .   Then  the  quantities 


(•■), 


J*)  -'    (\a('^\J'^> 


(rK     ,..   -_  /'\>^       ^"'^-(^> 


<<^'-=\-1^   ^7->     ('--■^>-) 


(7.9) 


satisfy 

(^^F,  (t)  =  ^"^F,  (tj  .  (7.10) 


■ka'  '       ka'  0' 


fr)  (r) 

We  find 


The  ^  '^F   (t)  may  be  evaluated  in  the  same  way  as  the    F  (t) 


^°)f,  (t)  =  b  ,  ^-^-ii) 

ka^ 


ka 
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.  b^d  I  ^)^^,\_,\v^_,%^(a)    ?  V^.3C2.2(«)]    >  (7.12) 

where  b  =   (e     +  l)      .      The  SLumnations  in   (7-12)   are  over  all  p  and  s  such 

that   I  pi ,     Is  L   and     k+p-s     are  all  less  than  k  fit   should  be  noted  that, 

since   H  is  finite,   the  problem  with  momentum  cut-off  admits  only  a  finite 

number  of  particles  per   system  <W)'/M-]      For  r  >  1,   we  find,   as  before,   that 

F      (t)    is  a  finite-degree  polynomial  in  b  and  in  the  C        (a)    of  order 
ka  n;p 

n    i    2r. 

(r) 
It  now  follows  from  (^-3)  that   F  (t)  has  a  bound  which  depends  on 

b  and  on  k  but  not  on  a,  t,  or  M.   Let  us  consider  the  complete  set  of  states 
s  such  that  W(s)/M  has  a  given  value.   Here  N(s)  =  \|^\^(s)  is  the  eigen- 
value of  the  total  number  of  particles.  The  boundedness  of    F  (t)  then 
implies  that  the  rth  derivatives  of  the  complete  set  of  matrix  elements 
/s'l  q.  (t)|  s^  have  a  mean-square  bound  which  is  independent  of  M  in  the  limit 
M  ->  oo  and  which  is  finite  if  W(s)/M  is  finite.   The  argument  is  the  same  as 
for  the  nonlinear  oscillator. 

7.3.   Integration  of  the  Heisenberg  Equations 

Let  us  now  consider  the  evaluation  of  the  nonlinear-oscillator  func- 
tions G  (t,t')  and  Q  (t,t')  by  a  finite-difference  solution  of  the  Heisenberg 
equations  of  motion.   We  seek  the  solution  over  some  given  time-domain 

(t^,  t   ).   Let  the  domain  be  divided  into  equal  intervals  At  and  let  (5.8)  and 
0^   max 

(5.11)  be  replaced  by  the  sets  of  difference  equations 
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and 

G  Jt     ,  ,t  )  =  G  ^(t  ,t  )  -  iAt  [eg  ^(t  ,t   )  +  M""Sm  ^(t  ,t   )1 , 

(l.lh) 
G      (t  ,t  )  =  6  „, 

where  t  =t+rAt   (r  =  0_,1^2^  ...)•   We  may  solve  these  equations  in 
terms  of  the  Schrodinger  operators  by  recursion,  substitute  the  results 
into  (5-19)  and  (5.20),  use  Wick's  theorem,  and  thereby  evaluate  the  quan- 
tities Q  (t  ,t  )  and  G  (t  ,t  ).   Equivalently,  we  may  substitute  into  (5.31) 

and  (5.32)  and  evaluate  the  quantities  J  (t  ,t  )  and  K  (t  ,t  ) .   The  re- 
^-^  -^   '  ^  a^  r-^  m       o  r  m 

suits  are  related  by  the  equations 

G  (t  ,,t  )  -  G  (t  ,t  )  =  A  tf-ieG  (t  ;t  )  +  J  (t  ,t   j]  , 
o^  r+1^  m'    O^  t'   m'  L    o^  r'   m'  O^  r^  m'J  ' 

Q  (t  ,  ,t  )  -  Q  (t  ,t  )  =  A  tf-ieQ  (t  ,t  )  +  K  (t  ,t   j]  , 
^o^  r+1^  m^    o^  v'   m'  \_         a^   r  m'    o^  v'   m'J  ' 


(7.15) 


which  are  the  finite-difference  forms  of  (5-27)  and  (5-29) • 

The  final  expressions  obtained  by  this  procedure  resemble  the  results 

of  the  iteration  expansion  carried  out  in  Sec.  6.   Each  contribution  to 

Q  (t  ,t  ) ,  or  to  the  other  functions,  consists  of  some  product  of  factors 
^o^  r'  m  •' 

Q   (tp^^t  )  and  Q    (tpi^t  )  multiplied  by  some  power  n  of  V  and  by  a  factor 

C   (o) .   As  in  the  case  of  the  iteration  expansion,  all  the  contributions 
n;p' 

containing  a  given  C  _  (o)  may  be  associated  with  the  pth  distinct  primi- 
tive linked  diagram  of  n  vertices.   There  are  two  differences,  however. 
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The  first  is  that  the  finite-difference  resiilts  involve  only  the  initial- 
value  Q   (t„,t_),  in  contrast  to  the  time-functions  Q   (t,t')  and 
G   (t,t')  which  appear  in  the  iteration  results.   This  is  because  we  are 
taking  finite-difference  approximations  to  (5.8)  and  (5-11)  rather  than 
to  the  integral  equations  (5-13)  and  (5.li|)  .  The  second  difference  is  that 
the  classes  of  diagrams  included  by  the  successive  stages  of  the  two  proce- 
dures are  very  different.  The  iteration  expansion  is  a  power-series  expan- 
sion in  V;  at  the  rth  stage^  it  yields  for  Q  (t,t')  an  approximation  which 

contains  only  the  C    (a)  of  order  n  ^  r.   In  contrast,  the  quantity 
n;p 

Q  (t  ,t  )  obtained  by  the  finite-difference  scheme  contains  some  C   (a) 
^a     r'   m  n;p 

as  high  as  n  =  2  (3""^  +  s"^  -  2) .  The  finite-difference  scheme  can  be  consi- 
dered a  particular  kind  of  consolidation,  re-ordering,  and  weighting  of 
the  iteration  (perturbation)  expansion. 

Let  us  now  consider  the  convergence  properties  of  the  finite-difference 
approximations  in  the  limit  A  t  ->  0,   The  exact  function  Q^(t,t')  may  be 
written  as  the  explicit  sum 

Q^(t,f)  =  Z-\^,exp[-aN(s)]<s|^^t')|s'><s'|^(t)|s>  .  (7-16) 

A  similar  expression  for  G  (t,t')  maybe  obtained  by  using  (5.IO)  .   The 
Heisenberg  equations  (5.8)  constitute  a  coupled  set  of  first-order  differen- 
tial equations  for  the  matrix  elements  which  appear  in  (7.I6) .   The  error 
in  the  matrix  elements  given  by  the  finite-difference  integration  scheme 
therefore  depends  on  the  magnitude  of  the  second  derivatives 
d^<^s'|  01  (t)|  s)  /dt^  over  the  domain  (t^,  ^^^y)  ■     These  derivatives,  however, 
are  not  boujided  for  all  s  and  s'.   Certain  matrix  elements  between  states  in 
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which  the  total  number  of  quanta  is  very  large  oscillate  with  extreme 
rapidity.   Wo  matter  how  small  A  t  may  he  taken,  there  will  be  matrix 
elements  that  are  poorly  approximated.   Convergence  of  the  finite-differ- 
ence approximations  for  Q  (t^f)  therefore  requires  that  the  weighting 
function  z"  exp[-aN(s)]  in  (7-l6)  suppress  the  contribution  of  matrix 
elements  with  infinitely  rapid  time-variation. 

We  found  in  Sec.  7.2  that  for  each  finite  r  the  quantities 


^^\it)    = 


a^\(t,t') 


satisfy  finite  bounds  which  are  independent  of  t,  a,    and  M.   We  note  also 
that  the  contribution  of  every  matrix  element  to  ^^'^F  (t)  is  real  and  non- 
negative.   It  follows  from  this  that  matrix  elements  with  infinite  time- 
derivatives  do  not  make  a  finite  contribution  to  Q  (t,t)  .   By  applying 
Schwarz's  inequality  to  (7.I6),  we  see  that  the  contribution  of  any  matrix 
element  to  Q  (t^f)  is  bounded  by  its  contributions  to  Q  (t,t)  and  Q  (t%t')' 
These  considerations  suggest  that  in  the  limit  A  t  ->  0  the  finite-difference 
approximations  to  Q  (t^f)  and  G  (t^f)  may  converge  to  give  the  exact  func- 
tions over  any  given  domain  (t  ,  t   ) . 

It  sho-uld  be  emphasized  that  we  have  not  given  a  proof  of  convergence. 
The  fact  that  matrix  elements  with  extremely  high  oscillation  frequencies 
make  a  negligible  contribution  to  the  ex^ct  Q  (t^f)  does  not  assure  that 
they  also  make  negligible  contributions  to  the  finite-difference  approxima- 
tions to  Q  (t,t').   The  successive  approximations  to  such  matrix  elements 
will  in  general  be  unstable,  and  can  greatly  exceed  the  exact  values  after 


63 


a  sufficient  number  of  integration  steps.   We  shall  not  attempt  to  resolve 
this  question  in  the  present  paper.  However^  we  conjecture  that  such 
instability  does  not  destroy  the  convergence  of  the  finite -difference 
approximations  to  Q  (t,t')-   The  reasoning  behind  the  conjecture  is  inti- 
mately connected  with  the  collective  nature  of  the  variables  qL,(t)  .   It  is 
clearest  for  the  true  problem,  and  we  shall  outline  the  argument  briefly 
for  this  case. 

The  Hamiltonian  for  the  true  problem  is  (5.I).   The  matrix  elements 
for  this  case  may  be  evaluated  immediately  in  the  representation  where  the 

[n]  "^[n] 
/s'l  qQ,(t)  i  s\   may  then  be  evaluated  by  transformation  to  the  representation 

with  the  a     a     diagonal.   If  one  carries  this  out  and  takes  M  very  large, 

the  results  strongly  suggest   that,  for  any  given  W,  those  /s'f  a  ,(t)|s)> 

which  have  significant  components  with  frequencies  »(e  +  W)  are  ones  for 

which  W(s)/(MN)  »  1.   However,  these  matrix  elements  are  very  strongly 

suppressed  by  the  weighting  factor  Z  exp[-aN(s)]  in  (7.I6)  when  M  is  very 

large.   This  makes  it  plausible  that  the  finite -difference  approximations 

to  Q  (t,t')  converge  as  A  t  ->  0  if  we  take  M  infinite.     Now  we  note  that 

the  C  _  (a)  are  independent  of  M  for  the  true  problem  and  consequently  the 

finite-difference  approximations  to  Q  (t,t')  are  independent  of  M.   This 

then  implies,  if  the  previous  argument  is  correct,  that  the  approximations 

to  Q  (t,t')  converge  for  any  M.   For  small  M,  however,  we  conjecture  that 

instabilities  arise  in  such  a  way  that  the  sum  over  states  s  and  s'  in 

(7-l6)  becomes,  in  effect,  a  sum  over  a  divergent  but  formally  correct 

series  as  we  take  A  t  ->  0. 
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The  argument  just  outlined  can  be  extended  to  the  stochastic  models  by 
introducing  a  'collection-of -collections'   consisting  of  M'  collections  each 
with  M  systems,  taking  collective  variables  in  the  collection-of -collect ions, 
and  considering  the  case  of  infinite  M'  ,   The  extension  can  also  be  made 
without  this  device.   The  central  part  of  the  argument  -  that  for  large  M 
only  states  with  W(s)/(mS)  »  1  give  rise  to  matrix  elements  having  signifi- 
cant components  with  frequencies  »  (e  +  VH)  -  appears  on  qualitative  grounds 
to  depend  only  on  the  collective  nature  of  the  q^(t)  and  to  be  as  valid  for 
the  models  as  for  the  true  problem.   A  point  of  consistency  which 
noted  here  is  that  the  explicit  expressions  for  the  successive  ^s'|  q^(t  )  |  s')  ^ 
obtained  by  recursive  solution  of  (7-13);  involve  only  matrix  elements 
between  intermediate  states  s"  which  satisfy  N(s")  ^  N(s)  -  1. 

It  now  remains  to  extend  our  considerations  to  the  many-body  problem. 
In  order  to  clarify  the  discussion,  let  the  weighting  operator  (5-35)  '^^ 
taken  in  the  particular  form 


W  =  bxp^-PqCHq  -  MqW)]  .  (7.17) 

If  V(x)  were  zero,  this  would  correspond  to  a  grand  canonical  ensemble  at 
temperature  P_  and  chemical  potential  M-  .  With  V(x)  nonzero,  (7-17)  does 
not  correspond  to  an  equilibrium  ensemble. 

We  may  now  write  Ql     (t,t')  as  the  explicit  sum 

X<s|C(*')l=='><s1\„(t)U>.  (7.18) 
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Here  s  and  s'  are  members  of  the  complete  set  of  joint -eigenstates  of  the 
number-operators  N   =  q  J'q   ,  the  N  (s)  are  the  eigenvalues  of  these 
operators,  and  W(s)  =  Z  N  (s)  is  the  eigenvalue  of  the  total-number- 
operator  N.   The  factor  Z  is  now  given  hy 


Z 


^(s)  -  Mq  N(s)1>  • 


v^  -poLVpV^^^  -^o5(^)J  •  (^-^9) 


If  the  term  -6^  2  e  N  ,(s)  were  ahsent  from  (7.I8),  our  discussion  of 
O  p7  p  p7 

the  convergence  of  the  finite-difference  approximations  for  the  nonlinear 
oscillator  would  be  immediately  applicable  to  the  many-body  problem  with 
momentum  cut-off.   This  follows  because  we  found  similar  mean-square  hounds 
on  matrix-element  derivatives  in  the  two  cases.   But  it  is  difficult  to  see 
how  the  presence  of  the  term  in  question  can  interfere  with  convergence. 
The  effect  to  to  weight  the  sum  against  states  in  which  there  is  strong 
initial  excitation  of  high  momenta  and,  therefore,  in  which  the  initial 
kinetic  energy  is  high.   Since  the  problem  is  conservative,  this  implies  a 
discrimination  for  all  t  against  certain  states  of  high  total  energy.   If 
we  consider  large  M,  then  the  weighting  operator  (T-I?)  effectively  confines 
the  sum  to  states  s  for  which  N(s)/(MN)  ^1,  where  N  is  determined  by  P^  and 
H  .   In  this  respect,  it  resembles  the  infinite-temperature  weighting  operator 
(7.8)  which  we  used  to  obtain  bounds  on  matric  elements.   But  further,  it 
chooses  from  among  the  complete  set  of  such  states  a  subset  which  is  weighted 
in  a  particular  fashion  against  high  total  energy  eigenvalues.   The  addition- 
al selection  may  reasonably  be  expected  to  increase  rather  than  decrease  the 
suppression  of  matrix  elements  with  extremely  high  frequencies  of  oscillation. 
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We  conclude  therefore  that  if  the  finite -difference  scheme  converges  for 
the  nonlinear  oscillator^  then  it  is  very  plausihle  that  the  finite-difference 
approximations  to  the  many-body  functions  a  (t^f)  and  G  (t,t')  also 
converge,  at  least  for  the  problem  with  momentum  cut-off. 

7.4.    Justification  of  the  Formal  Model  Equations 

Let  us  assiune,  on  the  basis  of  Sec.  7.3,  that  the  finite-difference 

approximations  to  G^^^(t,t')  and  Qj^^(t,t')  converge  for  all  M  as  A  t  ->  0,  and 

that  the  convergence  is  independent  of  M  for  large  M.   It  follows  directly 

that  the  total  contribution  to  the  unwanted  diagrams  vanishes  for  each  of 

the  stochastic  models  in  the  limit  M  ->  <».   The  unwanted  C   (a)  which  appear 

n;p 

in  each  approximation  are  of  finite  order  and  vanish  in  the  limit  M  -^0°. 
Thus  they  do  not  appear  in  the  final  functions  to  which  the  sequence  of 
finite-difference  approximations  converges. 

The  justification  of  our  formal  closed  model  equations  for  non- 
equilibrium  does  not  immediately  follow  from  the  vanishing  of  the  contribution 
of  the  unwanted  diagrams.   These  equations  were  obtained  in  Sec.  6  by  carrying 
out  formal  sums  of  the  primitive-diagram  iteration  expansions  for  a  (t,t') 
and  G  (t;,t')  with  only  the  wanted  diagrams  included.   However,  it  is  possible 
that  the  expansions  are  divergent  even  when  they  are  restricted  to  wanted 
diagrams.   In  this  case,  it  is  not  assured  that  the  formal  closed  model  equa- 
tions represent  a  unique  summation  of  the  expansions.   This  question  can  be 
resolved  by  considering  the  formal  model  equations  themselves  as  the  limits 
of  finite-difference  equations  rather  than  as  infinite  sums  of  perturbation 
terms.   The  justification  of  the  equations  proceeds  in  several  steps,  which 
we  shall  outline  in  terms  of  the  nonlinear  oscillator  example. 
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It  is  convenient  for  the  present  purpose  to  carry  out  a  finite- 
difference  solution  of  the  integral  equations  (5. 13)  and  (5-1^)  instead 
of  (5-8)  and  (5-ll).  Upon  substituting  the  results  into  (5-19),  (5-20), 
(5.31),  and  (5.32),  one  obtains  expressions  for  G  (t  ,t  ),    Q  (t  ,t  ), 

J   (t  ,t  ),  and  K  (t  ,t  )  which  may  be  evaluated  by  Wick's  theorem  in  terms 
a  r'  m^      a  r  m 

of  the  quantities  ^        i^^,>\)>    ^^  ^^*r ' '*m' ^  ^  ^^^*r"*m'^^  ^^^ 

G    (t  , ,t  ,).   These  expressions  are  more  compact  (and  also  more  acciirate 

^  r '   m' 

for  finite  At)  than  those  obtained  by  the  more  elementary  integration  scheme 

of  Sec.  7.3.   The  results  of  the  present  procedure  are  analogous  to  the 

primitive  expansions  of  Sec.  6.1.   They  contain  contributions  associated 

with  all  the  primitive  linked  diagrams.   Only  diagrams  of  finite  order 

contribute  to  G  (t  ,t  ) ;,  etc.  for  given  r  and  m. 

The  expressions  for  J  (t  ,t  )  and  K  (t  ,t   )  can  be  reformulated 
-^  arm      arm 

in  terms  of  contributions  associated  with  irreducible  diagrams  only.   As 
in  Sec.  6.2,  one  obtains  the  irreducible  forms  by  retaining  only  the  con- 
tributions associated  with  all  the  irreducible  diagrams,  and  in  them  replac- 
ing  each  factor  G^°^t^,  ,t^, ) ,  G^°)*(t^,  ,t^, ) ,  ^^^^  (*r  "^m"  ^ '  °^ 

Q  (t  ,;t  ,),       or        Q  (t  ,,t  ,),  where  o  is  the  collective  index 

which  labels  the  relevant  line.   The  irreducible -diagram  expressions  for 

J  (t  ,t  )  and  K  (t  ,t  )  may  be  verified,  as  in  Sec.  6.2,  by  substituting 
a^  r'   m      a^  r'  m 

for  each  factor  G  (t  ,,t  ,),  etc.,  its  finite-difference  expression  in  terms 

of  primitive  diagrams.   In  the  present  case,  the  expression  for  each 

G  (t   ,t   )  is  a  finite  sum  over  primitive -diagram  contributions  instead  of 
0^  r'   m' 

an  infinite  series.   Alternatively,  the  irreducible  expressions  may  be 
obtained  by  the  variational  method  mentioned  in  Sec.  6.2. 
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At  this  point  we  may  take  the  limit  M  ^  <»,  so  that  only  contributions 
associated  with  the  wanted  irreducihle  diagrams  survive.   In  the  case  of  the 
Hartree-Fock  and  random-coupling  models,  it  may  then  be  verified  that  if  the 
expressions  for  J  (t  ,t  )  and  K  (t  ,t  )  are  substituted  into  (7.I5)  the  re- 
sults are  just  the  finite-difference  forms  of  the  closed  integro-differential 
equations  for  these  models  given  in  Appendix  B.   In  the  limit  A  t  ->  0,  they 
become  identical  with  these  equations.   For  the  ladder  and  ring  models,  the 
further  step  remains  of  summing  the  (finite)  series  of  irreducible  diagrams 
which  contribute  to  each  quantity  K^("t^^"t^)  and  J^(t^^t^)  •   This  can  be  done 
by  introducing  vertex  functions  defined  by  difference  equations.   The  final 
sets  of  equations  thus  obtained  are  again  the  finite -difference  forms  of  the 
closed  integro-differential  equations  for  the  models  which  are  given  in 
Appendix  B. 

The  formal  closed  model  equations  may  thus  be  obtained  without  any 
use  of  the  perturbation  expansions  of  Sec.  6.   Their  justification  then 
depends  solely  upon  the  validity  of  our  assumption  that  the  finite-differ- 
ence integration  scheme  converges  in  the  limit  A  t  ->  0. 

7.5     Approach  to  Equilibrium 

Let  us  now  consider  the  extent  to  which  validity  of  the  formal  closed 
model  equations  for  nonequilibrium  implies  validity  of  the  equilibrium  model 
equations  of  Sec.  3.   On  the  basis  of  the  preceding  discussion,  we  shall 
suppose  throughout  the  present  Section  that  the  nonequilibrium  closed  equa- 
tions are  valid  over  any  given  domain  (t^^t   )  for  all  the  stochastic  models 
with  momentum  cut-off. 
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Consider  first  the  ladder  and  ring  models  when  the  conditions  stated 
in  Sec.  2  for  boundedness  of  the  potential  energy  per  particle  are  satisfied. 
It  is  easy  to  verify  that  the  imposition  of  a  momentixm  cut-off  does  not 
destroy  the  boundedness  properties.  The  interaction  Hamiltonian  for  the 
ladder  model  may  still  be  written  in  the  form  (2.3O),  which  is  positive- 
definite  if  V(x)  is  non-negative  everywhere.     When  all  the  V  are  non- 
negative^  the  lower  bound    "2  V(0)W,  implied  by  (2.33)  for  the  total  poten- 
tial energy  of  the  ring  model,  is  replaced  by  the  less  negative  bound 


-  t  s 


okp 


^   N,  ^,  where  the  sum  over  momenta  is  restricted  by  i  k I  ^  k 
k-p  ka'  J'  I   I     J, 


Ipl  S  k 
'   '     max 

Suppose  that  we  take  the  nonequilibriura  ensemble  in  the  form  (7.I7) . 

For  given  P^  and  ^^,   momenta  I  k  I  -^  k    will  have  negligible  initial  exci- 
^0      0'         '   '     max  '^  ^ 

tation  if  k    is  sufficiently  high.   Since  the  potential  energy  per  particle 
is  bounded  from  below  and  the  system  is  conservative,  it  then  is  plausible 
that  such  momenta  will  be  negligibly  excited  at  any  later  time.   We  anticipate 
that  the  behavior  of  G  (t,t')  and  a  (t,t')  will  be  independent  of  k  ____  for 


max 


the  nonequilibrium  fimctions  G  (t,t')  and  \Q,(*>'t')  should  depend  only  on 
t-t'  and  should  be  related  by  (5-^8)  to  the  temperature-domain  equilibrium 
propagators  for  some  P  and  U-.      (in  general,  we  will  have  p  7^  p_  and  H  ^  M._.) 

Since  the  unwanted  C  ,  (a)  do  not  contribute  to  the  non-equilibrium  functions, 

njp 

it  follows  that  they  should  not  contribute  to  the  equilibrium  propagators  in 
the  limit  M  ->oo.   The  closed  equations  for  the  model  propagators  are  then 
justified  provided,  in  addition,  that  the  formal  summations  of  wanted  diagrams 
described  in  Sec.  3  are  justified.   The  latter  supposition  is  made  plausible 
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by  the  discussion  of  Sec.  7.k,   which  validates  the  corresponding  summation  of 
wanted  diagrams  for  the  non-equilihrium  functions. 

The  arguments  just  presented  are  not  conclusive^  and  they  are  based 
in  part  on  ergodic  assumptions  which  are  very  difficult  to  investigate. 
However,  it  may  be  feasible  to  check  the  correspondence  between  the  equili- 
brium and  non-equilibirum  closed  model  equations  in  a  direct  analytical 
fashion.   If  both  sets  of  equations  are  valid,  and  if  our  ergodic  assump- 
tions (Sec.  5-5)  are  valid,  then  the  equations  for  S   (^  )  should  be  obtain- 
able by  analytic  continuation  from  the  none  qui  lib  rl  urn  equations  for  t  -^  -  <». 
We  have  not  attempted  this. 

Tlie  following  point  should  be  noted.   Some  equilibrium  values  of  p 
and  M  may  be  unreachable  by  any  choice  of  p  and  M-  in  (T-I?).   In  particizlar^ 
if  V(x)  is  purely  repulsive  the  initial  potential  energy  will  be  so  high  that 
very  low  temperatures  will  be  unreachable.   The  reason  is  that  (T-I?)  repre- 
sents zero  initial  two-body  correlations.   We  may  handle  this  situation  by 
taking  a  true  problem  in  which  the  system  of  interest  is  coupled  by  weak 
forces  to  a  reservoir  of  otherwise  free  particles  and  constructing  a  stochas- 
tic model  of  the  combined  system.   If  enough  particles  are  in  the  reservoir, 
and  the  coupling  is  weak  enough,  then  any  desired  equilibrium  temperature  may 
be  reached  by  evolution  of  an  ensemble  of  the  form  (T-I?) •   In  this  way  we 
may  justify  the  formal  closed  equilibrium  equations  for  the  propagators  of 
the  combined  system.   Finally,  we  may  let  the  coupling  to  the  reservoir  go 
to  zero  and  thereby  recover  the  closed  propagator  equations  for  the  system 
of  interest  in  isolation. 
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Now  let  us  consider  the  random  coupling  model,  in  which  there  is 

no  lower  bound  to  the  potential  energy  per  particle  in  the  limit  M  ->  <». 

Our  remarks  will  also  apply  to  the  ladder  and  ring  models  when  V(x)  does 

not  satisfy  the  conditions  which  give  bounds  on  the  eigenvalues  of  H. .   In 

-1 

the  absence  of  a  momentum  cut-off,  it  is  not  assured  that  these  models  will 
evolve  to  equilibrium  at  all.   It  is  possible  that  the  mean  potential  energy 
may  grow  negatively  infinite  and  the  mean  kinetic  energy  positively  infinite. 
With  the  momentum  cut-off,  there  is  a  ceiling  on  the  kinetic  energy,  and  in 
this  case  we  anticipate  that  an  equilibrium  will  be  achieved.   The  supposi- 
tion is  supported  by  the  existence  of  the  rigorous  equilibriiun  ensembles 
given  by  (7-8).   However,  such  an  equilibrium  may  differ  very  markedly  in 
its  properties  from  a  grand  canonical  ensemble.   In  particiilar,  it  may  be 
unstable  tmder  coupling  to  external  systems.   Obviously,  its  properties  need 

not  become  independent  of  k    as  k    ->«>.   These  considerations  lead  us  to 
max     max 

regard  the  formal  closed  equations  for  the  temperature -domain  propagators 
with  strong  suspicion  in  the  case  of  models  with  no  lower  bound  on  the 
potential  energy  per  particle.   The  classical  results  presented  in  I  perhaps 
cast  some  light  on  the  situation. 
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APEEEDIX  A.   GENERALIZED  MODELS  FOR  DISTINGUISHABLE  PARTICLES 

Let  the  true  Hamlltonian  for  a  system  of  N  distinguishable  particles 
interacting  through  the  pair  potential  V(x)  be 


H  =  J  £.  P.^  +  H.,  H.  =  ^  Z.  .'VCx.-  X.),  (A.l] 

2   1-^1     i'    i    2ij'i    j' 


as  in  I:  Sec.  2.1.   As  before,  L.   .'  means  that  i  =  j  is  to  be  omitted  in 
the  summation.   We  consider  a  collection  of  M  such  system  (M  odd)  with 
total  true  Hamiltonian 


H  =  J  Z  E  p.  ^,^  +  H.  ,    H.  =  i  Z  E.  .'V(x.  ^  ,  -  x.^   ),         (A. 2) 
-   2  n  i  i  [n]     -i'   -i   2  n  ij   '  i  [n]     J  [n]  ' 


where  x     and  p     are  the  position  and  momentum  of  the  1th  particle  of 
i  [n]     -^i  rn] 

the  nth  system.   Then  we  take  the  general  model  H.  in  the  form 


5i=i^nm^ij'^ijCnm]KM  "  ^  j  M  ^'  ^^'3^ 


V.  .  ^^^  (x)  =  M-lv(x)  Vxp  [-i2.(n-m)a/M]  ^.^..^  (A.k) 

[a  =  0,  +  1,  . . .  ,  +  I  (M-i)  ]  . 

To  recover  the  true  problem,  we  take  <^.     .     =   1  for  all  i,  j,   and  a, 
thereby  obtaining 

V.  .,  -,  (x)  =  8  V(x).  (A.  5) 

1 J  [nmj       nm 
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In  the  stochastic  models,  the  4-     -.rv  ^^-"ve  unit  modulus  hut  stochastically 
determined  phases.  In  this  case,  the  individual  systems  in  the  collection 
are  dynamically  intercoupled.   We  see  from  (AJ+)  that  ^^  ^  rnml  ^''^i  fnl  "■'^  i  fml  ^ 
depends  not  only  on  the  displacement  x    -  x     of  the  pair  of  particles 
in  space,  but  also  on  their  'displacement'  n-m  in  the  collection.  More- 
over, V.     (x)  is  different,  in  general,  for  each  pair  i, j. 

The  distinguishahle  particle  versions  of  the  generalized  ladder,  ring, 
random -coupling,  and  Hartree-Fock  models  are  given  by  the  following 
prescriptions: 

Ladder  Model 
Take 

-   exp 


i;j;a 


+i2rtaA.  ./m  ,   A  .  =  -A.  .  ,  (a. 6) 


For  each  pair  of  indices  i,j  fix  the  integer  A.  .  by  an  independent  random 
choice  in  the  interval  (O  :S  A  <  M) ,  subject  only  to  the  antisymmetry  con- 
straint in  (a. 6).   By  (A.U)  and  (2.10),  we  find 

V.  .,  Jx)  =  5  .  V(x),  (A.7) 

ijrnml       n-A.  .,m 

where  &   a      is  to  be  interpreted  according  to  the  cyclic  convention  (2.I5). 
n-^.  .,m 

Ring  Model 
Take 
h,y,a  -  '^   [-2.(A^.^  .  A.._^)/m]  ,   A^.^  =  -  A.  ._^.  (A.8) 

For  each  pair  of  indices  i,a  fix  the  integer  A.    by  an  independent  random 

1  ju 

choice  in  the  interval  (O  <A  <  m),  subject  only  to  the  antisymmetry  constraint 
in  (A.8). 


(A.9) 


Random-Coupling  Model 
Take 
</.     .        =  exp  -i2rt  A  .  M  \  , 

A      =  _  A      ,   A      =  -A 
i,j;a     j,i;c(         ±,j;a         ±,j;-a 

For  each  triad  of  indices  i; j^a,  fix  the  integer  A.  .   by  an  independent 
random  choice  in  the  interval  (O  —  A  <  m)^  subject  only  to  the  antisymmetry 
constraints  in  (A.9)- 

Hartree-Fock  Model 
Take 
5^.^..,  =  1,  ^,^^.^=0  (a  /O).  (A.IO) 

There  are  no  random  parameters.   By  (A.i|)  and  (2.10),  we  have 

V.  .^  ,  =  m"^V(x.  ,  ,  -  x.^  J.  (A. 11) 

ijrnm]         i[n]    j  [mj 


A  comparison  of  the  present  models  with  those  of  I:Sec.  2  shows  that 
' displacement -in-collection'  n-m,  and  the  associated  'Fourier'  modes  a, 
now  play  the  roles  in  constructing  the  randomized  potentials  that  formerly 
were  played  by  spatial  displacement  x.-  x.  and  the  Fourier  modes  k.   In 
the  present  models,  the  shape  of  the  interaction  between  any  pair  of  par- 
ticles is  always  that  of  V(x) ;  there  is  no  mutilation  of  the  potential  as 
in  I:Sec.  2.   However,  the  strength  of  this  interaction  can  vary  with 
n-m,  i,  and  j.   In  the  simplest  case,  the  Hartree-Fock  model,  we  see  from 
(a. 3)  and  (A.ll)  that  each  particle  simply  moves  in  the  average  field  of 
the  entire  collection  of  particles. 
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If  one  takes  a  grand  canonical  ensemble  of  collections,  the  present 
models  lead,  in  the  limit  M  -*  oo^  to  closed  expressions  for  the  classical 
Helmholtz  free  energy  whatever  may  be  the  value  of  N,  the  mean  number  of 
particles  per  system.   A  classical  non-equilibriiim  formalism,  analogous 
to  that  of  Sees.  5  and  6,  may  be  developed  for  these  models.   It  involves 

n-body  time -displaced  distribution  functions  (n  =  1,2,...)  and  also  Green's 
fxmctions  which  give  the  response  of  these  distributions  to  small 
perturbat i  ons . 

APEEWDDC  B.   NOMEQUILIBRIUM  MODEL  EQUATIONS 
We  list  here  the  final  expressions  for  J  (t,t')  and  K  (t,t')  for  the 
four  stochastic  models  in  the  limit  M  -^«>.   The  collective  indices  a,... 
are  omitted  since  there  is  no  dependence  on  these  indices  in  the  limit. 
The  equations  below,  taken  together  with  (5- ^3)  and  {^.kh),   form  complete 
sets  which  determine  the  evolution  of  G,  (t,t')^  "\(t,t');  and 
Ql     (t,t')  =  G  (t,t')  +  CI  (t,t').   (As  in  the  text,  the  upper  and  lower  signs 
of  a  double  sign  refer  to  fermions  and  bosons,  respectively. )  The  equations 
for  the  boson  random-coupling  model  (with  t^  =  -«>)  have  been  given  previously. 
The  equations  for  the  several  nonlinear  oscillator  models  may  be  obtained  from 
those  below  simply  by  omitting  all  the  momentum  indices  and  sums  over  momenta. 
All  integrals  in  the  equations  below  may  be  taken  from  t_  to  +^j  the  G  func- 
tions then  automatically  restrict  the  actual  ranges  of  integration  according 
to  the  defining  relation  Gj^(t,t')  =  0   (t  <  t'). 

Hartree-Fock  Model 

Jj^^(t,f)  =  -iZ^(VQ  +  V^__^)Np(t)G^(t,f)  (B.l) 

Kj^(t,f)  =  -i2^p(Vo  ^  V^_^)Np(t)Q^(t,t')  (B.2) 
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W  (t)  =  Q  (t,t)  (B.3) 

p     p 


the  mean  number  of  particles -per -system  with  momentum  p  at  time  t 


Random  Coupling  Model 

'Ss\J\-s^\-s^     I  [^G/(t,t")Q^(t.t")Qjt.f) 

+  Qp(t",t)G^(t,t")Q^(t,t") 

-  Q  (t",t)Q/(t,t")G  ft,t")'|G  (t",t')dt"  (B.k) 

p       r        b      J  -"^^ 


Vt,t>)  =  -i^p(Vo^v^,_p)Np(t)S,(t,f ) 

^  Vk-s(\-s^Vs^  I  [Gj(tSf)Q/(t",t)Q^(t,t")Q^(t,t") 


+  Q^(t",t')Gp  (t,t")Q^(t,t")Q^(t,t") 


(t",t' )Q  (t",t)G^(t,t")Q^(t,t") 


-  Q^^(t",t')Qp(t",t)Q^^(t,t")G^(t,t")jdt"  (B.5) 

In  these  equations  r  =  k  +  p-s. 


Ring  Model 

J^(t,t')   =   -iV^NG^Jt^t')  1  iZ^  |v^_^'(t.t^)Q^(t,.t^)G^(t^,f)dt^ 


-Z 
ps 


.ffi'v^.,'(t,t3)V,'(t,,t2)Q^(t,,tj)Q/(t3,t^) 

X  G^(t,tj_)Gj^(tj_,t'  jilt^at^dtj  (B.6) 
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X  G^(t,t^)]     dt^dt2dt3  (B.7) 


Here  N  =  2  N   (t)   is  the    (constant)  mean  ntmiber  of  particles  per   system,   and 


(t,t')    =  V^[s(t-f)   +   iE^,    _f[Gp,*    (t,f'   )   Qp,^^(t,f   ) 
-   S'^^"    '^^   %-+ci^^'*"   ^]     V^'(t'St')dt"j    . 


(B.8) 


It   satisfies 


V^'(t,f)  =  V_^'*(t,t').  (B.9) 

If  the  initial  momentiim  distribution  has  reflectional  symmetry,  V  '(t,t') 
is  real.   The  vertex  function  may  be  interpreted  in  terms  of  a  higher-order 
Green's  function,  of  a  kind  which  is  nonvanishing  only  for  nonlinear 
systems.   ¥e  have 

=  -'   ff\(*^^'\-s'(*''\' V^^'^'l^^r(^l''^2^ 

X  G  (t',t  )dt'dt  ' .  (B.IO) 

S       J        -L 
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Here  6f   (t )  is  the  Infinitesimal  perturbation  operator  introduced  in 
Sec.  5,   ^Rul'  ™sans  that  the  values  (^  =  a  and  X  =  a  are  omitted  from 
the  sum,  and,  again,  r  =  k  +  p  -  s.   Equation  (B.IO)  is  valid  only  for 
the  ring  model  and  for  M  — >  oo. 

Ladder  Model 

J^(t,f)  =  -iZ^  Jv^,^^^^'(t,t^)Q^(t^,t)G^(t^,f)dt^ 

^  ^ps  .QI%rs'^*^^3^ksr'*^\^^2V^'^'l^ 

X  Q^(t^,t2)Q^(t2,t2)Gj,(t^,t')dt^dt2dt^  (B.ll) 

+Z   \T\V    '(t,t,)V,    '*(t,,t^)  [g  *(t',tjQ  ^(t,  ,t) 
ps  .).).)  kprs  ^  ^  3'  kprs   '  1^  2'  *-  k  '   '  1'  p  '  1'  ' 

+  Q^(t^,t')G^*(t,t^)]  Q^(t2,t2)Q^(t^,t2)dt^dt2dt^      (B.12) 

Here  r  =  k  +  p  -  s,  and 

V    '(t,t')=V,    '(t,t')  +  V.    '(t,t').  (B.13) 

krps    '  kprs  kpsr 

The  vertex  function  V,    'is  determined  by 
kprs 

V,_...__'(t,f)  =  Vj._^&(t-t')  -  ±\,\_^,   |[+  G^,(t;t"  )Q^,(t,t"  ) 


kprs 


X  V^,^,^^'(t",t)  +  Q^,"'(t,t")G^,(t,t")V^,^,^^'(t",t)]dt". 

(B.lll) 
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It  satisfies 

V^    '  (t,t'  )  =  V  ,   '  (t,t'  ),  (B.15) 

kprs  pksr  "■    '      ^ '  \^--^^) 

a  fact  which  has  been  used  in  writing  (B.12).   In  analogy  to  the  case  of 
the  ring  model,  \   '  niay  be  interpreted  in  terms  of  a  higher-order  Green's 
function.   We  have 

^■V•W,^-^^%M^*4\a(^V[Sfp;(t,)6f^,(t2)&f3,(t3)]> 
=  -i  .{y\(t,t')V^^^^'(f,t2')G^Mt' ,t^)G^(t2',t2) 


y~Q^^{t^  ,t^)QX'^t^  ,  (B.16) 


where  ^o|a'  has  the  same  meaning  as  in  (B.IO).   Equation  (B.16)  is  valid 
only  for  the  ladder  model  and  for  M  ^  oo. 

If  the  vertex  functions  V  in  (B.6),  (B.?),  (B.II),  and  (B.12)  are 
expanded  into  infinite  series  by  iteration  solution  of  (B.8)  and  (B.lif), 
there  result  explicitly  all  the  contributions  associated  with  the  irreducible 
diagrams  that  survive  in  the  ring  and  ladder  models. 

The  nonequilibrium  equations  presented  above  for  the  several  stochastic 
models  should  reduce  to  an  equilibrium  description  in  the  limit  t^  ->  -<», 
subject  to  the  reservations  expressed  in  Sees.  5.5  and  7- 5-   In  this  case, 


the  equations  may  be  simplified  considerably  by  transformation  to  the  fre- 
quency domain  and  the  use  of  (5.51)-   The  results  possibly  may  prove  a 
useful  adjunct  to  the  model  equations  for  the  temperature-domain  propagators 
presented  in  I: Sec. 5*   In  some  applications,  the  equilibrium  quantity  of 
direct  interest  is  Q,  (w),  defined  by  (5'50).   In  principle,  Q,  (w)  may  be 
determined  from  the  temperature -domain  propagator  S  (^  )  by  analytical 
continuation,  but  in  practice  this  may  prove  very  difficult.   It  is  there- 
fore of  interest  to  have  equations  which  directly  determine  Q,  (w). 

In  this  connection,  it  should  be  noted  that  the  formalism  used  for  the 
temperature -domain  propagators  in  I,  and  in  Sec.  3  of  the  present  paper, 
does  not  appear  to  be  directly  applicable  to  propagators  in  the  real-time 
or  frequency  domain.   The  reason  is  that  the  individual  states  which  make 
up  the  unperturbed  (free-particle)  grand  canonical  ensemble  are  not,  in 
general,  eigenstates  of  the  adiabatic  S  matrix;  the  particles  can  scatter 
each  other  irreversibly  while  the  interaction  is  switched  on.   Consequent- 
ly, it  is  not  clear  that  the  perturbation  series  for  the  real-time  propa- 
gators can  be  expressed  in  terms  of  time-ordered  products  in  the  simple  way 
that  is  possible  for  the  temperature-domain  propagators. 
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FOOTNOTES 

1.  R.H.  Kraichnan,  New  York  University;,  Institute  of  Mathematical  Sciences, 
Div.  of  Electromagnetic  Research,  Report  HT-9  (196I). 

2.  R.H.  Kraichnan,  J.  Math.  Phys.  2,   12k   (196I). 

3.  The  Notation  I: (3-1)  denotes  equation  (3-1)  of  I.   Such  notation  will 
he  used  throughout  the  present  paper.   Sections  and  Figures  in  I  will  he  de- 
noted in  the  forms  I: Sec.  3  and  I: Fig.  2. 

h.        The  present  procedirre  should  not  be  confused  with  that  used  in  I: Sec.  6. 
There,  each  second-quantized  field  represented  a  single  particle  instead  of 
a  whole  many-body  system. 

5.  In  order  to  minimize  confusion  in  expressions  with  multiple  indices, 
we  shall  denote  individual-system  quantities  by  square-bracketed  Latin 
indices  and  collective-field  quantities  by  unbracketed  Greek  indices  through- 
out this  paper.   Momenta  and  'energies'  will  be  denoted  by  unbracketed  Latin 
indices,  as  in  I. 

6.  See  Sec.  3  of   reference  2  for  a  discussion  of  the  physical  significance 
of  this  similarity  in  a  classical  context. 

7.  Cf.  R.E.  Peierls,  Proc .  Roy.  Soc.  (London)  A2]J+,  II1.3  (1952). 

8.  For  the  true  problem  [Hamiltonian  (5.I)]  ,  the  exact  functions  are  also 
independent  of  a  for  N  ^  0.      For  the  fermion  case,  H.  vanishes  identically 


the  exact  results 


and  we  have  Q  (t,t')  =  Q   (t,t'),  etc.   For  the  boson  case,  we  easily  find 
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It  is  of  interest  to  note  from  these  expressions  that  if  Q  (t;,t')  is 
expanded  as  a  power  series  in  V,  the  radius  of  convergence  is  given  by 
|V|  =a/|t-t'  |. 

9.   The  methods  to  he  presented  in  Sees.  6  and  7  s^re   also  applicable  to 
(non-Gaussian)  normal  ensembles  which  describe  statistically  inhomogeneous 
systems.   Such  ensembles  are  specified  by  weighting  operators  of  the  form 


^  =  ^^  [  -^.pa^^^^'P^%.p^  V]  ' 


where  w(k;p)  is  a  suitable  function.   In  place  of  CI  (t;,t')  and  G   (t,t')^ 
one  must  deal  with  the  more  general  quantities 

Alternatively^  the  analysis  may  be  carried  out  directly  in  x  space  by  using 
the  fields  t  ,(x)  and  t   (x).   Such  a  treatment  is  illustrated  for  two  classi- 
cal field  problems  in  Sees.  10  and  11  of  reference  2.   It  is  also  possible 
to  work  with  non-normal  weighting  operators  which  describe^  for  example, 
non-zero  initial  two-body  correlations.   This  type  of  generalization  is 
briefly  discussed  in  Sec.  9  of  reference  2. 

10.  A  review  discussion  is  given  by  D.W.  Zubarev^  Usp.  Fiz.  Nau]-i  71,    71 
(i960).  [Translation:  Soviet  Physics  Uspekhi,  3,    320  (1960).]   Our  Gj^(t-t') 
differs  by  a  factor  i  from  that  commonly  defined. 
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11.  R.  Kubo,  J.  Phys.  Soc.  Japan,  12,  570  (1957). 

12.  P.C.  Martin  and  J.  Schwinger,  Phys.  Rev.  115,  13^2  (l959). 

13.  R.H.  Kraichnan,  Phys.  Rev.  112,  105lj-  (I958). 

Ik.      Cf.  C.  Bloch  and  C.  De  Dominicis,  Nuclear  Phys.  J,  ii59  (I958). 

15.  It  is  essential  here  that  the  Heisenberg  equations  of  motion  be  taken 
in  the  form  (5-8).   If,  instead,  the  equation 

d^dt  =  i(H^-  ^H) 

is  used  directly  in  the  iteration  procedure,  both  unlinked  and  linked  dia- 
grams arise,  just  as  in  the  Schrodinger  representation.  Thus  the  elimina- 
tion of  unlinked  diagrams  in  our  procedure  can  be  traced  to  the  use  of  the 
commutation  relations  (5-5)  in  obtaining  (5.8). 

16.  In  contrast  to  those  for  J  (t,t')  and  K  (t,t'),  the  primitive  expan- 
sions for  G  (t,t')  and  Q  (t,t')  have  factors  associated  with  both  external 
lines. 

17-  As  it  is  described  in  reference  2,  the  variational  method  is  applicable 
only  in  the  limit  M  -^  00,   However,  it  can  be  extended  to  finite  M  by  the 
formal  device  of  introducing  a  ' collection-of -collections ' ,  consisting  of 
M'  collections  each  with  M  systems,  and  considering  the  limit  M'  ^  00  . 

18.   A  nonequilibrium  linked-diagram  formalism  for  infinite  fermion  systems 
has  been  described  by  K.  Wishikawa  [j.  Phys.  Soc.  Japan,  I5,  78  (I96O) J  , 
who  uses  the  interaction  representation  instead  of  the  Heisenberg  represen- 
tation.  The  basic  quantities  are  taken  as  (1  (t,t')  and  Q,  (t,t'),  rather 
than  Q,  (t,t')  and  G  (t,t').   Equivalent  complications  in  the  construction  of 
primitive  and  irreducible  expansions  arise  in  that  formulation,  if  the 
analysis  is  carried  out  correctly  for  a  Gaussian  initial  ensemble  and 
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finite  t-t„.   [The  asymptotic  irreducible  expansion  expressed  by  Fig.  9  of 
the  cited  paper  Is  not  valid  for  finite  t-t^-J 

19.  We  have  not  proved  this  rigorously. 

20.  It  is  important  here  that  the  Heisenberg  equations  of  motion  used  in 
the  finite-difference  procedure  be  talven  in  the  form  (5.8)  and  not  in  the 
fonn 

dq^/dt   =   idq^  -  q^H). 

The  latter  form  yields  unlinked-  as  well  as  linked-diagram  contributions 
to  Q  (t,t'),  as  we  have  already  noted.   It  is  unsuitable  for  a  finite-dif- 
ference integration  procedure  because  the  eigenvalues  of  H  grow  with  M  so 
that,  for  large  M,   da  /dt  is  expressed  as  the  difference  of  two  operators 
each  with  large  matrix  elements.   In  this  situation,  it  could  not  be  expected 
that  the  convergence  of  a  finite -difference  procedure  would  be  independent 
of  M  as  M  -*oo.  An  equivalent  difficulty  arises  if  one  attempts  finite- 
difference  integration  in  the  Schrodinger  representation. 


Higher  V  are  admitted  where  they  connect  q  and  q  factors  admitted  by  the 
cut-off.   In  the  x-space  representation,  the  cut-off  leaves  V(x)  unaltered 
but  restricts  the  fields  t^(x)  and  t   (x)  to  Fourier  sums  which  do  not  con- 
tain the  excluded  momentum  modes. 


Fig.  1.   A  labeled  vertex. 
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Fig.  2.   (a)  The  diagram  associated  with  C   (a) ; 
(b)  the  diagram  associated  with  C  .(a). 


V  V  V  V 


< * < * < i <: • ^— 


Fig.  5.   Incomplete  linear  diagram  with  four  vertices. 


Q''\t,t) 


(a'^oc)^  (o^Wo^) 


Fig.  k.      Representation  of  the  primitive  contribution  to  K  (t,t') 
proportional  to  C   (a) . 

J-5-L 


Fig.  5-   Tl^e  diagram  as 


sociated  with  C   (a). 


G-'°'(t,t') 


Fig.  6.   Representation  of  the  primitive  contribution  to  J  (t,t') 
proportional  to  C,  ^,(a). 
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(a) 


(;)''\t",t) 


p-^^°'(t,t"; 


&'°'(t,t")         &'°VtV) 


(b) 

Fig.  7-   (a)  The  diagram  associated  with  C   (a) ; 

(t')^(c)^(d)  representations  of  the  primitive  contributions 
to  J  (t,t')  proportional  to  C   (a) . 

Note ;  See  following  page  for  Fig.  7,  parts  (c)  and  (d). 


Q^°^it'\t) 


(d) 


Fig.    1}   continued.      Parts    (c)    and   (d)  . 


NEW  YORK  UNIVERSITY 

INSTITUTE  OF  MATHEMATICAL  SCIENCE^ 

LIBRARY 

4  Washington  Place,  New  York  3,  N.  Y_ 


APR  6 

^^^?)ATE  DUE 

1 
1 

t^^sm 

1 

C.ORO 

10  Kraichnan 

!'T-10 

i\raichnan 

:iYi; 

HT-IO 
Kraichnan 


Stochastic  models  for  many- 
iTLE  body  systems...  ^t" 


-wa- 


gwT^I/rTSrifc 


BORROWERS 


N.  Y.  U.  Institute  of 
Mathematical  Sciences 

4  Washington  Place 
New  York  3,  N.  Y. 


